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ABSTRACT 

A state-of-the-art, tunable diode laser infrared heterodyne 
spectrometer has been designed and constructed for ground-based 
observations throughout the 8-12 micron atmospheric window. The 
instrument has been optimized for use with presently available tunable 
diode lasers, and has been designed as a flexible field system for use 
with large reflecting telescopes. The instrument was aligned and 
calibrated using laboratory and astronomical sources. 

Observations of SiO fundamental (v “ 1-0) and hot band (v = 2-1) 
absorption features have been made in sunspots near 8 microns using 
the spectrometer. These measurements constitute the first reported 
detection of solar SIO, and the first such measurement using coherent 
detection. The data permit an unajibiguous determination of the 
temperature-pressure relation in the upper layers of the umbral at- 
mosphere, and support the stinspot model suggested by Stellmacher aud 
Wiehr (Astron. and As trophys . , 432, 1970). 
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CHAPTER I 
INTRODUCTION 

1<1 Geneipal Statemgnt of the .Problem 

1.1.1 The Need for a General .Purpose,._Tunable .Diode Las9.f!_ 

H.sA9r,Q,dyjis-ip,gQ.ferjim9fe,gi: 

Infrared heterodyne speotrosoopy is a passive technique 

employing coherent optical detection for the study of spectral 

features in remote sources. It employs optical components such as 

mirrors and lenses normally associated with incoherent optics, but 

since it is a coherent process the principles of cm and mm wave 

receivers dictate the instrument design. Because of its eeherent 

nature, it offers the same advantages as the longer wavelength 

techniques, notably ultra-high resolving power and high frequency 

precision. Heterodyne spectroscopy can be considered complementary 

to, rather than in competition with, Fourier transform spectroscopy 

(FTS), which is widely used for high resolution IR measurements of 

remote sources. The FTS technique permits measurements of large 

5 

spectral intervals with frequency resolving powers v/Av < 10 . 
Infrared heterodyne spectroscopy covers a much smaller spectral 
interval in a single measurement, but can provide almost arbitrarily 
high frequency resolving power. This permits the recovery of fully 
resolved line shapes and reveals fine structure not accessible to 
the FTS technique. 

Since the advent of high speed infrared photomixers, heterodyne 
spectroscopy using CO gas lasers has become a powerful technique 
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for ’•.he of terrestrial atmospherlo spaoios and astrophysioal 

objects at high spatial and spectral resolution. Unlike coherent 
detection at radio frequencies, however, gas laser local oscillators 
only provide for operation at the fixed frequency of the laser 
transitions. Although a large number of transitions are available, 
the intermediate frequency coverage amounts to about 7 percent of 
the typical spacing between transitions, so these instruments only 
allow measurement over very small intervals throughout their 
overall operating wavelength range. Moreover, their operation is 
presently restricted to wavelengths between 8 and 12 microns. 

The substitution of tunable diode- lasers (TDL’s) for gas lasers 


aa local oacillators in heterodyne systems offers the premiss for 


continuous tunability, allowing operation over the entire nominal 
wavelength range of the diode laser. The development of a reliable 
TDL heterodyne spectrometer with the sensitivity of gas laser 
instruments has yet to be realized, however, because of the many 
problems associated with the development of tunable lead-salt diode 
lasers. Instead of lasing in a single gaussian mode, the output 
power of these devices is typically divided among a mixture of 
high-order modes, and the output exhibits a variety of noise 
problems related to multi-mode emission. Also, because of their 
tunability, the output frequency is sensitive to mechanical 
vibrations, small temperature fluctuations, and small variations in 
the diode laser injection current normally used to fine tune the 
frequency. 

Aside from these difficulties, tunable diode laser heterodyne 
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spectroscopy has been successfully demonstrated In measurements of 
laboratory gases, atmospheric constituents and continue of the moon 
and Mars, Meanwhile, improvements nave been made in the performance 
and reliability of tunable lead-suit diode lasers, particularly with 
regard to power output, 

In view of the success achieved with gas laser heterodyne 
systems, and the demonstrated potofftial of the diode laser 
heterodyne technique, there is an obvious demand for a general 
purpose ground based tunable diode laser heterodyne spectrometer, 
for both astrophysioal observations and for measurements of 
terrestrial atmospheric constituents. 

jMga3»r.gfflgnt^a£,SiOJ.n.,Sun g ,gfl. t , a 

One very interesting class of observational problems that is 
ideally suited for TDL infrared heterodyne spectroscopy is the 
detection of molecular absorption in sunspots. Since temperatures 
in sunspots are lower than that of the surrounding photosphere, the 
environment is suitable for the formation of a variety of molecular 
species, some of which have already been observed. The sunspot 
environment is of itself interesting since it is not well understood 
and the roles of magnetic fields, laaoroturbulence and heating from 
the overlying chromosphere in determining the environment have not 
been well established. In addition, the sunspot environment 
constitutes an excellent high-temperature source for the study of 
species not easily created in the laboratory. 

In order to properly analyse absorption lines observed in the 
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sunspot spectrum, it is essential to Know the vertical 
temperature-pressure distribution throughout the region where line 
formation takes place. Unfortunately, wide disagreement exists 
between empirical one-dimensional sunspot temperature/ pressure 
models that have appearea in the literature. The disagreement is 
largest in the upper layers of the umbral atmosphere where the 
formation of most atomic and molecular lines occur. 

For the deeper layers, models employ continuum observations 

at both visible and infrared wavelengths. One difficulty that 

exists with sunspot continuum observations is contamination by 

terrestrial lines and the unfc®): ainty in the transmittance of the 

terrestrial atmosphere. Since atmospheric transmission depends on 

local conditions at the time of measurement, these effects are 

removed by recording spectra from the spot and photosphere in close 

time coincidence and ratioing to obtain the relative sunspot 

intensity. Also, at visible wavelengths the ratio of photospheric 

to sunspot surface intensities is large (I , /I > 10), so that 

ph urn 

observations are contaminated by scattered light from the 
surrounding photosphere. Each continuum measurement must include 
some method of correcting for this effect. Attempts to model the 
upper layers of the umbral atmosphere have been based on 
observations of atomic lines at visible wavelengths where the 
sunspot spectrum is very complicated. This makes analyses of 
individual line shapes exceedingly difficult because of blending 
with adjacent lines. 

Most of the above-mentioned problems are not as serious at 


qoauty 
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infrared wavelengths. The ratio I .71 deoroasea at longer 
wavelengths and the continuum speotrura shows discrete, isolated, 
well defined absorption features, few continuum measurements have 
been made at wavelengths longer than one micron, however, and 
observations of infrared absorption features are almost 
non-existent. All infrared line and continuum measurements are thus 
of groat importance in the study of the sunspot environment. 

The fundamental (v = 1-0) and hot band ( v c 2-1) transitions 

of silicon monoxide (SiO) near 8 mioronr are of special interest. 

These transitions have been predicted in the sunspot speotrum but no 

detection has bean reported in the literature. Also, SiO is a 

sensitive thermometer at the upper layers since it is rapidly 

dissociated at sunspot temperatures. By comparing the observed 

equivalent widths against calculated values based on assumed 

temperature - pressure profiles it should be possible to determine 

the temperatures at small optical depths and provide boundary 

* 

conditions for improved models. The tunable diode laser heterodyne 
technique is ideally suited to these observations for several 
reasons. The 8 micron SiO bands occur at the frequencies of present 
oomraeroially available diode lasers. Also, since it is a coherent 
detection technique, the field of view is diffraction limited, with 
a diameter of about 3 arseo at 8 microns using a 1 mater telescope 
aperture. This beam size is smaller than most urabral diameters so 
contamination- from surrounding photospherio light is minimal , at 
least during periods of good seeing, Since the technique is capable 
of subdoppler resolving power, it permits line broadening and 
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residual intensity measurements, in addition to the usual equivalent 
wiath information. Finally, the high surface brightness of sunspot 
umbrae (T st 4000K) makes SiO observations a logical first choice 
among the list of astrophysically interesting problems for the TDL 
heterodyne technique. 

1.2 f rogress in IR Heterodyne Detection Using 

Dig,<a,Q,,Laagn,..LogaI,.QsQj.l3,aV?rs 

Prior to the development of milliwatt semiconductor tunable 

diode lasers and wideband HgCdTe photodetectors, sensitivity 

limitations in diode laser heterodyne detection work required the 

use of strong signal sources. Success in early heterodyne studies 

was achieved by reporting observed self-beating effects in a single 

source or by detecting optical beats from two IS local oscillators. 

Hinkley et al. (1968) observed beating effects between a current 

tunable PbSnTe diode laser and the P(20) transition of a COg gas 

laser. The beat note at approximately 1250 MHz was observed using a 

* 

copper doped germanium (Ge:Cu) detector, with 40mW of CO^ laser 
power and approximately 100 uW dingle mode power from the diode 
laser. The instantaneous linewidth of the diode laser emission was 
found not to exceed some tens of KHz, although frequency modulation 
of the diode output was observed cue to a 120 Hz ripple in the 
current supply. These authors also observed heterodyning with the 
P(18) COg transition at IF frequencies below 3.5 GHz, with somewhat 
poorer amplitude and frequency stability. 

Heterodyne measurements of the Lorentzian spectral distribution 
of output power from a liquid helium-cooled PbSnTe tunable diode 
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laser were perforaed by Hinkley ana freed (1969) in an attempt to 
verify theoretical predictions for the dioae laser line shape. 
Because of the intrinsioally small cavity dimensions and low output 
power levels of TDL's, diode laser linewidths are significantly 
greater than those of gas lasers operating in the same frequency 
range, so that the gas laser line can be oonsiderec a delta function 
by oompa\rison. By combining both diode and CO^ gas laser radiation 
at a liquid helium cooled Ge;Cu detector and current-tuning the 
diode output frequency across the extremely narrow CO^ laser line 
(Av < 100 Hz), the resulting beat note was detected and the entire 
diode laser line displayed on a spectrum analyzer. The 
experimentally determined linewidths were found to agree with 
theoretical predictions based on estimates of the diode "oold" - 
cavity bandwidth and the degree of population inversion. The 
authors were also able to demonstrate the inverse power variation of 
the diode laser linewidth, by increasing the diode current and 
adjusting the heat sink temperature to maintain a constant beat 
frequency. 

The first successfully operating infrared heterodyne 
spectrometer using a tunable semiconductor diode laser local 
oscillator was constructed by Mamma et al. (1975) and used to make 
laboratory heterodyne measurements of line profiles in the 
band (100-000). These authors also demonstrated heterodyne 
detection of themal blaokbody emission from the moon and Mars. The 
spectrometer employed a PbSe semi-conductor diode laser, a HgCdTe 
photodiode as a photcmixer and a spectral line receiver consisting 
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of 8 contiguous 25 Mriz filters to provide a coverage of 200 MHz at 

the IF. At the operating wavelength of 8,6 microns, this 

8 

corresponds to a resolving power of about 1.A X 10 , NgO heterocyne 
absorption measurements were performed using a cell filled with N^O 
at 10 torr pressure placed in the beam of a 1300 K blaolcbody 
continuian source. The line width was about 170 MHz, which was 
fully resolved in the bank of filters. 

Blaokbody heterodyne measurements of the moon and Mars yielded 
measured signal to noise ratios of 6.7 and 5.4, respectively, for a 
single 25 MHz filter and roughly 8 minutes integration time. Much 
of the signal degradation vras the result of insufficient local 
oscillator power. The maximum single mode power obtained in this 
experiment was >^100 nW, considerably below that required for shot 
noise limited operation (i.e., the point where the local oscillator 
shot noise dominates all other noise sources at the IF) . From the 
measured photodiode quantum efficiency and diode laser single mode 
power, the overall transmission of atmosphere and system optics was 
estimated to lie between 10 and 20 percent. 

Frerking and. Muehlnar (1977) developed a diode laser heterodyne 
spectrometer and obtained laboratory and terrestrial atmospheric 0^ 
absorption spectra in a 1 cm”"* region near 1011 cm”’^ This 
instrument employed a PbSnSe laser diode and HgCdTe photodiode mixer 
both mounted, together with the optical elements, in compact fashion 
on a common dewar operating at liquid nitrogen temperature (77K) . A 
single IF amplifier and filter at the mixer output with a cutoff of 
100 MHz provided a single IF channel with a width of 100 MHz, 
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ylQldins an IR resoivia? power of about 1.5 X 10"^ . Spectral scans 
of the band (001 000) of atmospherio ozone in solar absorption 

were obtaineo by slowly ourrent-tuning the diode between 1010.9 and 
1011,8 om"^ while recording the synchronously detected heterodyne 
signal on an X-Y plotter, lixoellent agreement was observed between 
the heterodyne spectra and direct absorption laboratory spectra 
using an ozone absorption cell. The results also agreed nicely with 
synthetic 0^ spectra in this region computed by the same authors, 
and adjusted to match the instrumenc frequency resolution., the 
measured SNR was found to be about 25 for a 100 MHz IF bandwidth and 
1 sec integration time. This is roughly a factor of ildOO smaller 
than the ideal value. The degradation was attributed to optics 
transmission losses, insufficient local oscillator power, amplifier 
noise and spurious IF noise. 

The first use of a closed cycle cooler for continuous diode 
laser cooling and temperature control in a heterodyne system was 
reported by Ku and Spears (1977). These authors constructed a 
heterodyne radiometer using a tunable 10.6 ym PbSnSe diode laser ad 
a local oscillator. The closed cycle cooler held the diode near 30K 
for the duration of the experiments. A heterodyne signal was 
obtained by combining the local oscillator beam and that from a 
1173K blackbody source at a HgCdTe photomixer. The heterodyne 
detection electronics consisted of a 600 MHz low noise amplifier at 
the mixer, a spectrum analyzer, crystal detector and look-in 
amplifier as part of the synchronous detection process. The 
instrument was tuned by varying both cooler temperature and the 
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diode current , Initial etuaies of noise at tne It’ suggested that 

self-beating effects in the laser output were the cause of extreme 

noise in the heterodyne signal over portions of the TDL tuning 

range. For quiet operation, signal- to-noise ratios as high as 150 

were observed for a blaokbody temperature of 1173K, a bandwidth B of 

500 MHz and an integration time t of 0.4 sec. In addition one of 

the laser modes was examined by heterodyning against a COg gas 

laser, and an instantaneous linewidth of 100 kHz was measured, in 

qualitative agreement with earlier results of Hinkley and Freed 

(1969). The instrument was also used to measure several absorption 

lines of ethylene near 942 cm in a gas absorption cell, The lines 

were first located in direct absorption and then observed in the 

heterodyne mode. ITie measured line widths were in good agreement 

with those expected for an IF resolution of »/> 600 MHz, 

A hybrid spectroscopic technique employing both direct 

absorption and heterodyne measurement options available using diode 

* 

lasers has been employed by several workers for the study of 
molecular gases in the laboratory. By monitoring the beat frequency 
between a PbSnSe diode laser and various COg gas laser lines, 
Worchesky and ooworkers (1978) were able to establish precise line 
center positions of several DgO absorption lines in near coincidence 
with COg gas laser transitions. Both the diode and COg laser beams 
were focussed onto a HgCdTe photomixer while monitoring the 
heterodyne signal to establish a frequency calibration scale. The 
remaining part of the diode laser radiation was passed through the 
DgO absorption ceil and onto a second detector where the 
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tranamittanoe waa aimultaneoualy monihoreci, A to*;al of nine linea 

in the v_ fundamental bands of 0.0 were measured with an estimated 
< 2 

error of ±0.005 om“\ Instabilities in the diode laser output 
frequency were the largest contributor to this uncertainty. 

Hillman et al. (1973) used a similar approach to study tne 
2sQ(1,1) transition in the band of near 9^9.5 cm"\ The 

spectrometer employed a PbSSe tunable diode laser mounted in a 
closed cycle cooler and emitting near the 9.4 am and 10.6 pm COg 
laser bands, Examination of the TDL output on a spectrum analyzer 
revealed an anomalously large mode width of »/’ 15 MHz resulting from 
radiation losses in the diode stripe, diode injection current noise 
and cooler instabilities. The mode center position stabilized to 
within ±2 HHs over several minutes of tise after the diode 
temperature had reached equilibrium. A small additional mode was 
also seen in the observed heterodyne spectrum about 60 MHz from the 
main mode, The presence of this type of satellite mode in the TDL 
output was thought to be the cause for regions of extreme excess 
noise of the type observed in heterodyne studies by Ku and Spears 
(1977). After an initial study of instri,fe;ent performance, a number 
of scans of the 2sQ(l,1) line were made using a 30-cm cell of NK^ at 
1-Torr pressure placed in one of the divided TDL beams. Aside from 
small departures from linearity in the TDL tuning rate, data from 12 
scans yielded a determination of the line ceiiter position to within 
± 0.00007 cm”^ or ± 2.1 MHz. Spears and Freed (1973) examined the 
use of IR - microwave varacting for continuous heterocyne coverage 
between adjacent CO^ transitions near 10.6 microns. This technique 
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also permits accurate viavelenssth calibration of a tunable diode 
laser using a COg laser as a reference. Usiag a tunable microwave 
local oscillator pump and a HgCdTe detector, it is possible to 
detect optical beat frequencies much larger than the single 
conversion IF bandwidth of the detector. By injecting the output of 
a microwave pump with frequency f^ into the mixer, the optical beat 
frequency sources is down converted to a final 

frequency f = nf^ - (fg - f^) which lies well below the nominal 
detector bandwidth. Spears and Freed (1973) measured a signal-to- 
noise exceeding 45 dB at 21.3 GHz with a 10 icHz noise bandwidth 
using two COg lasers operating at different transitions and 
dellivering ^ 0,5 mW each to the detector. The signal-ho--nQiae 
ratio was limited mainly by the available microwave local oscillator 
power. Subsequent 3rd and 4th harmonic mixing studies (n = 3,4) 
resulted in the detection of beat frequencies between 40.3 GHz and 
60.7 GHz. 

* 

The most recent, and probably the most aur.cessful attempt to 
employ TDL heterodyne spectroscopy for the study of atmospheric 
constituents was made during 1979 by Harward and Hoell (1980), The 
authors observed atmospheric NH^, HgO, 0^, COg and HNO^ absorption 
between 9 and 11 microns using a TDL radiometer and a small 
heliostat located at NASA Langely Research Center, Hampton, VA. The 
system employed a TDL local oscillator in a closed cycle cooler, the 
output of which was collimated using an f/1 Ge lens. Mixing df the 
LO and signal beams was accomplished with a HgCdTe photoqJixer. 

Coarse frequency calibration was performed using a monochrociiator , 
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and prior to eaoft observation the final local oaoiiiator position 
was sat usincs reference gas line positions and a solid Ge efealon to 
provide a frequency scale. One interesting operational feature of 
the instrument was the use of a scanning etalon to select and track 
a single longitudinal diode laser mode. 

The atmospheric spectra were recorded by varying the TDL 
current to sweep a single 100 MHz IF filter in frequency. Signal- to- 
noise ratios of 200 to 300 were obtained during these measurements. 
This study established the capability of the TDL heterodyne 
technique for terrestrial atmospheric studies between 9 and 11 
microns. 

1.3 Tunable Di,gila.JL,aag.tiJ?aLS.I.a.gtne,nt. 

Two major classes of tunable semi-conductor diode lasers exist 
today - the first based on the III-V binary compound GaAs, where 
lasing typically occurs in the 0.6-5 micron spectral range. The 
Reman numerals refer to the number of valence electrons for the 
constituent atoms and defines the role - donor or acceptor - of the 
constituent. The state-of-the-art has advanced rapidly for this 
class of device largely because of their potential for optical 
communication applications. Present GaAs devices are very reliable 
and exhibit multi-milliwatt power output, hign efficiency, beam 
uniformity and room temperature oparatit,^. 

This section vd.ll be concerned exclusively vd.th the second 
major diode laser class which is based on lead-salt, or IV-VI 
compounds. Interest in these devices stems from their potential for 
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wide tunability, their operation in the Important 3-30 raioron ranje 
(which includes the vibration-rotation bands of many molecular 

-4 

species), and their extremely hi?h spectral purity (less than 10 

-1 -1 

cm at V >/' 1000 cm ) , 

Butler et al. (1964) first reported lasins in the Pb-salt 

semi-conductor family using a PbTe (lead-telluride) diode laser 

operating at 6.6 microns. Lasing throughout the 2.7 to 33 micron 

region was later achieved after the development of ternary compounds 

employing basic binary constituents. 

The nominal wavelength operating range of a Pb-salt diode laser 

is obtained during fabrication by the proper choice of material, 

which in turn determines the energy bandgap. Figure 1 shows the 

variation of energy bandgap E for various ternary Pb-salt compounds 

S 

and the nominal operating wavelengths of these devices. Because of 

thj ffect known as band inversion, E approaches zero as x 

E 

approaches 0.4 in the compound Pb^ ^Sn^Te. In principle, this 
allows fabrication of very long wavelength devices. Several methods 
also exist for fine tuning the operating frequency of these devices. 
The output frequency can be smoothly varied by changing the applied 
pressure, the surrounding magnetic field, the injection current 
flowing through the device the mount temperature. The latter two 
techniques have proven the most popular since they are easily 
performed under field operating conditions. 

Early tunable diode lasers (figure 2a) were broad area 
homojunction devices employing a P and N type material and a narrow 
active region. This region acts both as a barrier to confine 
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minori*’.y oarpiora (the bandgap changes afc the sides of the junction) 
and as a radiation waveguide (a change in the refractive index 
occurs at the boundaries) , The front and back; surfaces are cleaved 
parallel and sometimes polished to define a Fabry-Perot cavity, The 
sides are roughened ,to minimize reflections and the generation of 
unwanted high-order modes, 

The operating spectral range of TDL's attained by temperature 
and bias current tuning can be specified, but the operating 
frequency range is much more restrictive because of mode hopping 
effects. The resonant mode frequencies in a device of cavity length 
I, and index of refraction n are given by tne Fabry-Perot condition 

Mil 

Ves “ ®/2nh cm (1,1) 

The cavity is usually many wavelengths long so "ra" is typically a 

large integer. When the diode is fine tuned by varying the 

* 

injection current, changes in the Joule heating of the material 
alter the temperature dependent index of refraction and hence the 
optical length. The effect of thermal expansion on the diode tuning 
is negligible by comparison. By subtracting Equation 1,1 from 
itself for adjacent "m", one obtains an expression for tne frequency 
separation Av of adjacent modes, 



dv 


Av s [2nL (1 + 


n 


(1 . 2 ) 
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Typical dimensions for lead-salt diode lasers are 400 pm length and 

200 (iffl in width and height, with typical output mode separations of 

'/'2 cm”^. Changing the device temperature also tunes the energy 

bandgap E in the material, typically several times more rapidly 
E 

than V for Pb-salt TDL’s, This differehtial motion changes the 
res 

preferred laser mode when the diode is current tuned, sometimes 

causing a mode “hop", and very often results in large frequency gaps 

over which no lasing occurs. In general the gain envelope, as 

determined by E at the operating temperature, has a complicated 
S 

shape and is wide enough to support several modes. This results in 
multi-mode output (Fig. 3) and necessitates provisions for external 
mode selection in many applications. The intent in fsibricating a 
TDL is usually the generation of a single, powerful longitudinal 
mode having TE polarization. Most performance problems are related 
to poor optical confinement in the laser cavity which creates higher 
order unwanted modes. These modes sliare the available power between 
them, which causes an irregular output power pattern and limits the 
tunability. Early ternary devices had high threshold currents, 
lased only at the lowest temperatures and operated with parasitic 
bounce modes involving all sides of the laser chip. 

A substantial improvement in parasitic mode suppression was 
realized using a stripe geometry for improved carrier confinement 
and mode guiding. In Figure 2b, the n layer formation is restricted 
to a narrow stripe on the p-type substrate using an SiO^ diffusion 
mask, and a metallized contact is finally applied. The mode guiding 
cavity along the junction direction is bounded by the ttripe edge. 
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and aidebounce modes are effectively attenuated by losses in the 
bulk material . The first use of a stripe geome.try was reported by 
Ralston et al. (1973) » fabricated using Pb. ooSn. ,„Te. Sidebounce 
modes were effectively suppressed and the emission spectra of the 
stripe consisted of regularly spaced frequency modes corresponding 
to a well defined longitudinal cavity. 

In addition to bounce-mode suppression, the performance of a 
TDL also depends on the degree of confinement of the injected 
carriers in the direction of the carrier drift, perpendicular to the 
junction. This is accomplished by the formation of heterostructures 
or heteroj unctions. These confine the carriers and the emitted 
light with epitaxially grown regions of wider energy gap. This 
decreases the lasing threshold current and reduces the resistive 
heating, permitting operation at higher amoient temperatures. This 
result is very significant since TDL operation at or above liquid 
nitrogen temperature (77K) simplifies the diode cooling apparatus. 
Single heterojunction devices formed by vacuum deposition (Walpole 
et al., 1973) and by compositional interdiffusion (CID) and 
liquid-phase-epitaxy (LPE) (Linden et al., 1977) exhibited low 
threshold currents and high operating temperatures, showing 
improvements in internal quantum efficiency. These features 
inspired the development of double heterostructure lasers for high 
temperature operation. Groves et al. (1974) achieved high power 
(^1.2mW) single mode cw performance in a stripe geometry diode laser 
operating at 77K. This is shown in Figure 2c. The active region 
was formed by LPE growth of n-type PbSnTe through 50 ym wide stripe 
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openings in a MgFg maf/k. Diffusion into the PbTe substrate was 
prevented by doping the substrate with T1, and the top border of the 
v/>10 )m wide active region was defined by the interface with a final 
PbTe layer. 

The fabrication of heterostructures is hampered by defects that 
form in the heterojunctions and in the surrounding material. These 
result from the difference in lattice parameters between adjacent 
layers. Fonstad et al. (198O) have successfully dealt with this 
mismatch problem using multi-layer ternary compounds having nearly 
identical lattice parameters. The threshold current (at 8OK) of 
these devices is a factor of >/'20 smaller than that for the 
PbSnTe/PbTe structure, and pulsed operation has been achieved to 
^166K. 

As discussed earlier, there is a tendency for a "mode-hop" to 
occur since the peak optical gain in the cavity tunes more rapidly 
with temperature than do the Fabry-Perot modes. This occurs even 
for a well defined cavity where bounce modes are suppressed and only 
the fundamental transverse mode exists. By replacing the 
Fabry-Perot end mirrors with a distributed feedback (DFB) cavity it 
is possible to greatly suppress this mode-hopping tendency (Walpole 
et al., 1976, 1977; Hsieh and Fonstad, 1979). This arrangement 
(Figure 2d) employs a weak periodic variation in the refractive 
index of the active region along the direction of wave propagation, 
with one end of the cavity beveled to suppress the generation of 
Fabry-Perot modes. In a DFB device the longitudinal modes are 
determined by the periodicity along the cavity and the single 





22 

allowed mode will be the one closest to both the Bragg frequency and 
the center of the optical gain envelope. Hsiah and Fonstad (1979) 
achieved single mode emission continuously tunaole over '/'20 cm”^ 
near 12.8 microns. Thia level of performance was attained without 
incorporating some of the previously discussed innovations (e.g., 
lattice matched layers) , so future efforts in DFB development should 
yield extremely wide continuous tunability. 

For heterodyne and other high signal to noise applications of 
TDL's, the level of excess noise present in the output is an 
Important consideration. Several broaa classes of diode excess 
noise related to multi-mode emission exist. These extend in 
frequency from several Hz to over 1 GHz and have bersn studied by 
various authors (Ku and Spears, 1977; Jennings and KilLman, 1977; 

Eng et al., 1979; Harward and Hoell, 1979). The noise sources 
include broadband noise, noise generated by mechanical vibrations 
associated with closed cycle cooler operation, uncontrolled feedback 

I* 

from optics external to the laser cavity, and the effects of mode 
competition within the diode. The various classes of high frequency 
TDL noise as documented by Harward and Sidney (1980) are shown in 
Figure 4. Most restrictive in heterodyne applications are the 
spiking resonance symptom (b) which seems to be associated with the 
heterodyne beating of spurious frequency modes, and broadband excess 
noise (e) covering the entire IF spectrum. Because of these 
effects, only a small fraction of the nominal diode tuning range is 
suitable for heterodyne work. 
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Rgure 4, Categories of diode laser excess noise (Harward and Sidney, 1980). 
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1.4 fjis.ylgjAs,.aunaflg|^-Ms, a M n .8,ig, n !; .3 ,,. ^ 

Most previous identifications of raoLeouiar absorption features 

in the sunspot spectrum have led to statements concerning the 

temperatures and pressures in sunspots. At the very least, attempts 

to observe molecules in spots imply some assumption about the 

ambient temperature. The converse is also largely true. Modelers 

of the temperature and pressure distribution in sunspots employ data 

on the size of molecular features along with continuum intensity 

measurements at several wavelengths. These data provide boundary 

conditions which, together with the rules for model stellar 

atmosphere construction, allow the computation of the vertical 

temperature and pressure structure, A review of sunspot models, 

then, must also include a description of sunspot observations, 

Wohl (1971) made one of the first synoptic studies of umbral 

features. The observations were made in 3 large umbrae during June 

and July of 1969 and covered the region from 400 to tiOOnm. The 

> 

instrument used was a photoelectric spectograph with a frequency 
resolving power of >^300, 000. Previous published molecular 
concentration studies were used by Wohl as a ledger in searcning for 
spectral lines. Because of photospheric scattered light and the 
effects of line blending, many of the lines could not be confirmed. 
Following a search for some 15,000 features, the following species 
were found to be present based on agreement with their expected 
positions and strengths; C^, CH, CN, CaH, i^H, TiO, CoH, HF, NH, 
and NiH. Tentative detections of the following species were also 
reported: AlH, BO, BeH, AIF, BF, BH, BaO, CuH, FeO, LaO, MgO, ND, 
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SoO, VO and ZrO. Because of a laoK of gooa molecular constants, it 
was not possible to accurately predict the frequencies and 
intensities of selected lines. This limited the number of 
identified molecules. It was also not possible to observe many 

weaker features since the instrument signal- to-noise ratio limited 

-4 

the detectable equivalent width to > 1 X 10 nm. 

Two of the early one-dimensional sunspot models, employing both 

infrared continuum measurements and observations of medium and 

strong atomic lines in tha visible spectrum, were published by 

Henoux (1968, 1969). Observations of umbral continuum intensities 

in a large spot (r >10 arc sec) were made at seven wavelengths 

between 0.36 and 2.2 microns using a 9 meter spectrograph and 50 cm 

telescope. Defining t as the 0.5 micron continuum optical 

0.5 

depth, Henoux' s (1968) model for t > 0.3 was computed from these 

0.5 

measurements by inverting the radiative transfer equation. 

Additional atomic line observations measured between 4?0 and 640nm 

* 

at a dispersion of 0.02 nm/mm were reported by Henoux in 1969. 

-4 

These enabled him to extend the model to roughly -t. _ j'lO . The 

0.5 

temperature structure for was adjusted to fit wings of 

observed Na D lines which change very rapidly with temperature. 

After obtaining a plausible temperature profile, the measured 
equivalent widths for unblended lines of Fell, CrII, Till, Fel, CrI 
and Til were used to check for aiireement with the model by varying 
the microturbulence parameter. A micro turbulent velocity of 1.5 
km/ sec yielded good agreement with observations. All of the line 
intensity measurements could then be corrected for stray light by 
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oomputing expected visible continuum intensities, basec on the more 
reliable infrared values. 

Only a few detections of sunspot molecules have oeen made at 

infrared wavelengths. (Hall and ooworkers at Kitt Peak), The 4 

micron SiO sunspot spectrum was examined by Hall in 1969 but the 

results remained unpublished (Hall, 1981). Hail and Noyes (1969) 

did report measurements of seve.'.'al unbral HF and comparison CO lines 

in both umbra and photosphere near 2,4 microns. These observations 

were made at the McMath solar telescope using a vacuum spectrograph 

with a resolution of 150,000* Although the intent was to estimate 

the solar fluorine abundance, obtaining this quantity from the 

observed absorption spectrum requires a fairly accurate temperature 

model, since hydrogen fluoride dissociates rapidly at sunspot 

temperatures. Effective temperatures for each spot were obtained by 

measuring the ratio of spot to adjacent photosphere intensity and 

comparing with absolute photospheric continuum estimates and model 

* 

limb darkening relations. The temperatures deduced for two large 
spots using this approach was 4000 1 200K and for the small spot , 

4400 ± 250K. The authors generated sunspot models consistent vdth 
the observed effective temperatures using the SAO stellar 
atmospheres and line synthesis program. The observed CO lines have 
widely different lower level energies so they were used as selection 
criteria for the choice of the best temperature model. After 
adopting a suitable model, the fluorine abundance was varied to 
obtain agreement with the measured equivalent widths. Hall and 
Noyes deduced a fluorine abundance (A^) of ■ 4.56 ± 0.33 on a 
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scale where ■ 12. 

Sfcellmaoher and Wlehr (1970) tested several previously 
published models , including Henoux (1969) using photographic and 
photoelectric measurements of 4 magnetically non-split lines, which 
were mostly free of blends. The authors then constructed a new 
model for the upper Ictyers, Lines of fa 1406.5402 nra, 543 **1527 nm, 
Fb’*’ 1722.4464 nm and Ti 1571.3895 nra were measured and continuum 
intensity measurements were made in line-free windows at 1461.5 nm, 
542.7 nm, 615.7 nra, 703.5 nm and 778.7nra for comparison purposes. 
Parasitic light corrections for those measurements were estimated 
using the ratio of observed lines of Fe 1614.925 nm, which is 
weakened in spots, to V 1615.015» which becomes stronger. 

Stellraaoher and Wiehr computed wavelength and center to limb 
vafiationo of the continuum as well as profiles of 3 of the measured 
non-split lines using published sunspot models and assuming LIE. 

Best agreement with both the continuum and line measurements 
occurred for the Henoux model, which was adopted for log ^ > 0. 
The authors also cited the previous range of sunspot microturbulent 
velocity determinations ( 1.1 - 2.9 km/ sec, Bray and Loughhead 
(1964); 1.8 km/sec photosphere, 1.2 xm/sec umbra, Bruckner (1965)), 
and quote a measured value of < 1.0 km/ sec by fitting the model to 
' the Fe 1406.54 line. A value of J’0.5 Km/sec, which is considerably 
smaller than moat previous solar microturbulent velocity 
determinations, was deduced from comparison with the Ti and Fe 
1543.45 nm features. 

An attempt was made by Webber (1971) to construct a one 
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dlmanaional sunspot model using a large body of visible absorption 
features obtained from 4th and gth order spectrograms taken at Mt, 
Wilson. Since the effects of scattered light varied greatly during 
the measurements the anphasis was on a rotational analysis of the 
lines within each spectrogram. Equivalent widths of several hundred 
observed lines of TIO, MgH and CaH were measured as part of Weboer's 
rotational analysis. The CaH analysis was oomplioated by the 
presence of excessive line blending and terrestrial atmospheric 
absorption. Different rotational temperatures were deduced for each 
species, allowing an assessment of the vertical temperature 
structure when Zwaans (1968) sunspot model was used to approximate 
the depth of iine formation for each moleouie, The final 
temperature profile was obtained by applying a small perturbation to 
the sunspot model of Henoux (1968). Some disagreement seems to 
exist between Webber's proposed model and the effective temperatures 
obtained from the very data used to construct the model, Also, the 
use of two separate previously published models in the analysis 
(Zwaan, 1968; Henoux, 1968) suggests a possible lack of self- 
consistency in the results, 

Kneer (1972) published a sunspot model based on both 
photoelectric measurements of the visible sunspot continuum at two 
wavelengths, and high resolution photographic spectra of 3 non-split 
Fel lines (X512.37, 543*45 and 557*6lnm) in 3 large umbral regions. 
These observations were corrected for stray light by extrapolating 
aureole intensities measured photoelectrically. After adjusting a 
hydrostatic equilibrium model to force the computed lines and 
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continuum intenaitiec to agree with the obaervations, the author 
produced a model in aatiafactory agreement with the reauita of 
Stellmaoher and Wiehr (1970). Since the winga of the obaorved Fel 
linea are only weakly preaaure dependent, the model preaauroa were 
uncertain by a factor of «/'2, ao amall departurea from hydroatatic 
equilibrium could not be ruled out, Thia atudy auggeated a 
mioroturbulent velocity of < 1 to^/aeo, in rough agreement with 
Stellmaoher and Wiehr (1970), anil an eatimate of -^1 kra/aec for 
umbral maoroturbulenoe. 

Some of the aupporting raeaaureraenta used by the previously 
cited authors were summarized by Zwaan (1974) in an attempt toj (i) 
explain discrepancies between publi? sd models and (ii) construct a 
seffli-empirioal continuum model atmosphere which reproduces measured 
continuum intensities between 0.48 and 2.3 microns. This study 
emphasized IR continuum measurements between 1,5 and 4 microns. 
Supporting observations were taken from Henoux (1968), Hall (1970), 
Coupiac and Koutchmy (1972), and Maltby v\nd Staveland (1971). The 
overlapping meaaurements of Henoux (1968) and Hall (1970) were 
basically in agreement, lending support to the model in the deepest 
layers, near the 1.7 micron H" opacity minimum. Additional visible 
and near IR photometric measurements were used to determine the run 
of temperature in the shallower layers. The model was fit to the 
visible data using a wavelength dependent opacity enhancement 
factor, since line blanketing is a dominant source of error at 
visible wavelengths. 

Figure 5 shows Zwaanfe observational summary and the variation 
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Temperature vs monochromatic optical depth at X - 0.5 um for some umbtal model atmo. 
spheres. The positions of unity optical depth at various wavelengths are shown by vertical ban The 

wavelengths arc indicated in ;tm. ' 


Figure 5, Summary of sunspot continuum measurements and published 
umbral models prior to I97I4 ( 2 waan, I97I*). 
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of surface Brightness with X, oal45ulated using several previously 

cited models, figure gb shows the models tneraselves and the 

discrepancies, which are minimized near log t_ _ * -0.3 to 0.0. 

u*p 

Zwaan's result is consistent witn radiative equilibrium for t. _ < 

0 . 0 

1.7 and departs to smaller lapse rates for t. _ > 2. This model is 

w.b 

considered reliable between 0.5 < t < 2 but has some limitations. 

W • p 

The run of temperature is given out to ^ .01, but the data 

# 

give no check on the model for t. . < O.S. The author also poi.its 

0 • P 

out the existence of a temperature inversion somewhere at small 
optical depths because of lateral radiative influx from the 
photosphere and heating from the chromosphere above the sunspot 
fluxtube. The modeig of Henoux (1969) a.nd Stcllmacher and niehr 
(1970), on the other hand, have some experimental justification at 
smaller optical depths, since they are based on fairly strong 
atomic line observations. The elevated temperatures in these models 
probably show the effects of heating from adjacent regions. 

Zwaan's (1974) model temperatures were also criticised as being 
too high for t > 1.5, since the Harvard Smithsonian fieference 
Atmosphere (Gingerich et al., 1971) used to convert the measured 
umbral intensities to brightness temperatures, was also considered 
to be too high (Vernazza et al., 1974). In response to this, Zwaan 
revised the deeper layer temperature structure, lowering the lapse 
rate to just above the adiabatic value (Zwaan, 1975). It should be 
pointed out that this change does not alter the validity of the IR 
observations used by Zwaan; only the conversion of the observations 
to umbra! temperatures are affected. 
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Hall and lloyes (1972) observed a number of rotation - vibration 
lines of in the sunspot spectrum as part of an effort to map 

the 2400 - 3000 cm"^ region. The measurements were made at the Kitt 
Peak McMath solar telescope using an IR spectrometer with a 
resolving power v/Av > 10® and a signal to noise ratio > 100. The 
frequency uncertainly for the measurements was «/> .005 cm"\ The 
authors detected 10, v=1-0 transitions and 4, v=2-1 transitions for 
H®®C1. Equivalent width measurements were made for all except 4 of 
the v=1”0 transitions, which were blended with adjacent lines. The 
results yielded a solar chlorine abundance log^^ N(C1) = 5.5 ± 0.4 
for both isotopes, on a scale where Ios^q ^(HIsIP. 

Gaur and coworkers published a number of separate studies 
investigating the presence of selected molecules in the sunspot 
environment. By performing dissociation equilibrium calculations 

s. 

using Zwaans (1965) sunspot model and using published estimates for 
elemental abundances and dissociation coefficients, Gaur et al . 
(1973) produced a search list containing the expected concentrations 
of a number of molecules in sunspots. The average formation depth 
for molecules included in these equilibrium calculations was also 
deterained. The molecules considered include those already observed 
plus many not yet observed at the time of publication (1973)* Table 
1 is a reproduction of the search list published by Gaur et al. The 
list does not include CoH and NiH discovered by Wohl (l97Di since 
Gaur did not have dissociation constants for these molecules. 

Gaur ( 1976 ) also computed the equivalent widths for selected 
Av=1, R branch transitions of NO near 5.3 microns in sunspots. The 
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coaiputationa were done for J = 27 using *:he sunspoH models of 

Stellmacher and Wiehr (1970) and Henoux ( 1969 ). The dependence of 
the computed equivalent widths (W) on solar disk position was also 
examined by calculating W at solar zenith angles (p) of 1, 0.7, 0.5 
and 0 . 3 . For each model, W was found to increase with p but 

systematically larger W's were computed using Henouxi model. The 

-4 

computed values ranged from 10 < W < 30 X 10 nm. 

A similar study (Qaur at al.,1978 ) was performed for SiO by 

calculating equivalent widths in the 8 micron sunspot spectrum for 

v=1-0 and v=2-1 P and R branch transitions near J s43. These 

max 

calculations were done with Zwaans (1974) sunspot model using 

published values for the SiO dissociation energy and equilibrium 

constant. The ^ formation depth fo.' all of the computed lines 

-4 

was about 0.063 and the equivalent widths ranged from 40-95 X 10 nm 
as was varied from 8.1 to 8.33 eV. The computed widths were also 
not very sensitive to the lower level vibrational quantum number. 

The line formation depth and temperature sensitivity for 
dissociation, however, show the usefulness of SiO measurements as a 
means of determining the temperature structure at higher layers. 

This study by Gaur et al. formed the incentive for the SiO 
observations discussed in Chapter V. 
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CHAPTER II 

THEORY OF INFRARED HETERODYNE DETECTION 

2>i Iha.,B a aA.flJi g lgr.gjtong..0.at§>Q,t^ ££g.aai8.a 

The theory of infrared heterodyne detection is well understood 

and has been reviewed by a number of authors (Blaney, 1975; Cohen, 

1975; Abbas et al., 1976; Betz, 1977). The intent here is to 

provide a brief but concise review of these fundamentals as they 

apply to an operating heterodyne system. 

Figure 6 shows the basic IR heterodyne process in schematic 

form employing an infrared local oscillator and a monochromatic IR 

signal source. Coherent detection occurs at a high spaed 

pho tods tec tor/mixer where the output current varies as the square of 

the total electric field incident on the detector. If (x,y,t) 

and S- (x,y,t) are the complex local oscillator and signal electric 

field contributions at a detector of active area A , the 

> 

intermediate frequency (IF) current from the detector will be 

ijj (t) s C E*(x,y,t) • E(x,y,t) dxdy (2.1) 

^D 

where 

E(x,y,t) = E^Q(x,y) exp(-iu»j^^t) + Eg(x,y) exp(«iWgt) 
and C is a constant. The surface integral in Equation 2.1 reveals 
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Figure 6. 5?he Ijaaic heterodyioe detection process. 
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the jjnportanoe of proper alignment in ‘an operating system so that 
phase fronts from the local oscillator and signal beams are properly 
matched across the detector surface. Also, the dot product shows 
that only a portion of the randomly polarized signal combines with 
the laser local oscillator to produce coherent detection. For a 
randomly polarized signal source, this means a factor of two 
reduction in the available signal power (Section 2.3). After 
separating the time dependence in Equation 2.1 and taking the real 
part of the result, one obtains 


yt) = 0 A„ 




'] 


+ 2C cosCui^^Q - 0 . 3 ) t |Ej_^q|1E31 


f 




( 2 . 2 ) 




In Equation 2.2, e, -'^nd g„ are the unit vectors for E and E . 

LO 5 uU u 

The constant C can be eliminated by examining the DC term in 
Equation 2.2^ which represents the detector's incoherent response to 
the total power incident on the active area. For a deceotor with 
constant quantum efficiency n, the optical power to IF-current 
conversion constant is ne/hv, e being the electronic charge. This 
implies 


|E 


LO' 


“l-q 

hv Ap C 


1/2 


(2.3a) 


and 


38 


OHIGIIMAl PAGE fS 
OF POOR QUALITY 


- no Pj, ~| 1/2 

hv A. C } (2.3b) 

L a 


so that in terns of the LO and source powers, Equation 2.2 can be 
written 


[' 


i„(t) = ^ I Plo * ^ 

hv 


(2.4) 


with 


S(t) 


1 

X 


cos (“lo-^sX' 


(e 


LO 


Cg) dxdy 




The first two terns in Equation 2.4 represent the photocurrent i^^^ 

arising from the total IR power on the detector, coherent or 

otherwise. For a randomly polarized signal source, the integral of 

* 

®L 0 “ I //27 and the intermediate frequency signal term 

ijp(t) in the expression for the detector current ip(b) becomes 

i (t1 = (P p ^ Qoa (u) -u_H (2.5) 

3-Ip ^ LO s' 

with a mean square value of 

< lj/> = (ne/hv)2 (2.6) 

after the integration over time. 
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2.2 Charaoberlsblos of th e Source Radiation 

In reality, the source relation for nearly all natural IR 
sources (e.g. the terrestrial atmosphere, planets, nebulae) will 
also be incoherent in n i' ure and distributed over a broad band of 
wavelengths. The usual problem entails observations of absorption 
or emission features against a blackbody continuum at a temperature 
T, in which case the background flux can be regarded as constant 
over the IF bandwidth and equal to the Planck function B;^(T) at that 
temperature. Even when observations of non-thermal sources are 
involved, the frequency width of the emitting feature is usually 
greater than a typical IF resolution element so it is use'/ul to 
define a brightness temperature T_ for any source. The intermediate 

D 

signal frequency in Equation 2.5 is then the center frequency 

Lw o 

of the IF bandwidth element of width B__, and the signal power P_ is 

Xr o 

the power per unit bandwidth integrated over B^.^.. For small B^^, 
the single sideband power P (watts) becomes 

w 



The terms A and n in Equation 2.7 refer to the aperture area and 
field of view of the collecting optics, and are required in addition 

•ip —1 —I 

to bandwidth B__ to convert B (T) in W m” Hz“ ster to received 
IF v' ' 
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power. For a coherent receiver, the field of view (fl) of the 
oolleoting optica ia juat the diffraction limited field of view for 
a collecting aperture of area A, The antenna theorem (Siegman, 
1966) impoaea an additional conatraint relating A, 0 and the 
wavelength 1, given by, 


A n X‘ 


( 2 . 8 ) 


The product Ah ia called the etendue, and is a constant of the 
entire system when all components are properly aligned and matched. 
This quantity, applied to the active photodetector area A^, sets the 
precision required in the alignment to guarantee that the signal and 
local oscillatar beams are properly phase matched across the 
detector. Equating with Ah, Equation 2.7 becomes 


P3 = 




(,2.9) 


so that the mean square IF current ia 


<i3./> = 2 

hv 


V B, 


( 2 . 10 ) 


Since noise sources must be taken into account before the 
instrument performance can be estimated, of great interest in a 
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heterodyne syetern is the signal- to-noise ratio (SNR) . The dominant 
noise source arises from statistical fluctuations in the arrival of 
photons, primarily from the local oscillator, and is called shot 
noise or sometimes quantuum noise (Oliver, 1965). The mean square 
shot noise current is determined by the detector DC photoourrent 
(ijjj,) and IF bandwidth and is given by 

<ig^>^ = 2 e igg s '’s' ®IF 


Normally >> Pg, so the shot noise current is local oscillator 
dominated and is simply called LO shot noise. 

If sufficient local oscillator power is available so that the 
LO shot noise dominates all other noise sources, the signal-to-noise 
ratio will be maximized and the expression for the predeteotion 
signal-to-noise ratio (SNR) takes a particularly simple form 


SNR 



( 2 . 12 ) 


This case is referred to as shot noise limited operation and shows 
that the signal-to-noise ratio becomes independent of local 
oscillator power at sufficiently large P . The dependence of 

uU 

signal-to-noise on IF bandwidth and integration time appears when 

the IF output is square-law detected and time integrated. For an 

integration time of t seconds and an IF bandwidth in Hz, the 

1 /2 

signal to noise ratio improves as (Bt) so that the post-detection 


A2 


value l3 given by 



for shot noise limited single sideband operation. 


(2.13) 


In oases where the available LO power is l,lmited, it is 

necessary to consider noise sources of thermal origin which occur as 

a result of the physical temperatures of both the mixer and the 

preamplifier. For a matched system with characteristic impedance 

2 

R^, the mean square noise current contributed from these items 

is related to the sum of their equivalent noise temperatures by 


2 

*> 

A 




4k 


^^A V ®IF 


(2.14) 


Most mixers currently used for IR heterodyne work are mercury 

cadmium telluride (HgCdTe) and operate at liquid nitrogen 

temperatures (>/’77K) in order to reduce their contribution to the 

total noise current i.. 

A 

When the amplifier and mixer noise contributions are taken into 
account, the post detection signal-tp-noise ratio has the more 
general form, 


SNR 


<lj/> . <l/> 


(2.15) 
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Using the above expressions for <isH^> and 

multiplying by 2 for double sideband operation > and making the 
approximation Pj^q P^, the signal-to-noise ratio becomes ■ 



( 2 . 16 ) 


which reduces to the shot noise limited case (Equation 2.13) when 
P, . > (hv/ne^)(2k/fl- H t. + T..) . In reality, the cest 
anti-reflection coated HgCdTe detectors have quantum effioier.-jies of 
about 0.8 at low frequencies, decreasing to 0.2 - 0.3 above 1 GHz. 
The dynamic impedance for these devices is usually greater than 
100 Q. The amplifier noUe temperature la determined largely by 
the first amplifier stage after the mixer. Wide-band low noise uhf 
amplifiers have measured noise figures typically less than 2.0 db 
throughout the 0.1 - 2,5 GHz IF range ofansidered here. This 
corresponds to a noise temperature of < I 70 K. Using these 
representative values, Equation 2.16 shows that laser local 
oscillator powers of 1Q0-200UW are required for near shot-noise 
limited operation at 8 microns. Lead-salt tunable diode lasers 
exhibit comparable single mode powers, and under optimum conditions 
can provide operation near the shot noise limit. 

In practice, the performance of an operating heterodyne system 
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la alao degraded by loaaea in opfcioal oomponenta, iraporfeob raatohlng 
of the algnal and looal oaolllator beams and as a result of the 
aynohronoua deteotlon process used in separating the signal frca the 
noisy background. The loasea associated vdth tunable diode lasers 
are the most severe because of their diverse noise properties and 
their marginal power output. 

Equation 2.13 nay be modified to include all sources of system 
performance degradation (including n) by inserting the term delta 
(A) in the denominator, viz., 


SNR s 


2 (B^p. 

A [exp (hv/kT)-1 ] 


(2.17) 


The term A is an effioienby parameter introduced by Abbas et al, 

(1976) which represents the degradation in SNR from all sources of 

system noise and loss of- heterodyne signal. By definition, A =1 for 

ideal shot-noise limited, single sideband operation with no 

* 

polarization or chopping loss, and increases as losses are 
introduced. In practice, a factor of 4 in degradation always exists 
since only one polarization component of the source radiation, that 
aligned with the LO field, is used while heterodyning, and the 
signal is chopped with a SO-SO duty cycle, A total a as small as 
7 has been achieved by the GSFC CO^ laser heterodyne system now 
operating at the Kitt Peak McMath solar telescope (Kostiuk, 1981). 
For diode laser systems, however, the A is usually higher, as 
discussed in Chapter IV. 

An alternate way of expressing the heterodyne sensitivity is in 
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teraa of noise equivalent power (NEP) . This is defined as the 
source power per unit frequency which produces unity signal-to-noise 
at the output. For a system operating at frequency v, the noise 
equivalent power is 


NEP a 


10 hv A 


W Hz 


-1 


(Bjp T) 


1/2 


(2.13) 


The factor 10""^ in Equation 2.18 converts cgs laiits to watts. 
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CHAPTEfl III 

* 

DEVELOPMENT OF THE SPECTROMETER 

3.1 Reoulremenhs for a High SenglUvlfcv. 8-1g .Micron 

TDL Heterodyne SDeohromef.er 

In order to produce atafce~of»t;he-arh performance from a fcunable 
diode laser heterodyne instrument, the overall design must 
incorporate features to reduce or eliminate the problems discussed 
in Chapter I. These problems are mostly related to, but not limited 
to, the multitude of noise characteristics associated with TDL's and 
the problems of environment and frequency control for these devices. 
To minimize these difficulties and produce a flexible field system 
for general asfcrophysical use, the instrument should have the 
following features; 

(i) A Qr;fogeinic£lly cooled, vibrationally isolated 
environment for the TDL local oscillator, capable of continuous 
long-term cooling. This minimizes device degradation resulting from 
repr'«-ed thermal cyclings to room temperature, which tends to occur 
with He dewar systems. 

(ii) High precision temperature control (i. 0.3 mK) of the TDL 
for periods of one hour or more, as required for extended 
integration times on weak sources. 

(iii) Provisions in the optical layout for matching the local 
oscillator beam to a wide range of signal beam focal lengths. This 
feature is essential if the instrument is to be moved to more than 


one field location. 
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(iv) Inclusion of a mode selector grating to eliminate all but 
one of the Fabry-Perot modes emitted simultaneously by the TDL, 

(v) Incorporation of diffraction limited reflecting optics 
wherever possible for improved instrtment throughput and to avoid 
the need for frequent realignment as the operating wavelength is 
changed. This problem is encountered with refractive optical 
elements , 

(vi) A dual-beam synchronous detection scheme capable of using 
any combination of remote source, remote reference (sky), local 
source and local reference in the signal and reference beams. 
Synchronoukt detection is essential for the extremely low heterodyne 
signal levels, and the choice of signal and reference must be 
flexible to accomodate a wide variety of observational approaches. 

(vii) Multiplex and time integration capability to extend the 
weak signal detection limit. 

3 . 2 EnvlrQnment..and-. »r e.q.uencv_..C.ontrol.jrnr._bhe..DlQde_._ 

The local oscillator presently installed in the instrument and 
used for nearly all of the observational work was a compositionaily 
inter-diffused (CID), PbSe tunable diode laser (TDL #8300-14), 
manufactured by Spectra Physics, Inc., Laser Analytics Div. 

Referring to Figure 1, the composition of this diode lies at the 
right on the vertical axis defined by the binary material Pbse. The 
output from this laser is generally multi-mode with a total power 
that varies between lOOpW and 1 mW, depending on operating 
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frequency. The power in a single mode, as required for heterodyne 
work, is generally leas than 300pw, however. The diode operates 
from about 1170 om"^ at a diode current of 0.20 A and mount 
temperature of 12 K, to nearly 1260 cm"^ at 2.0 A and 60 K. 

Proper mounting of the device poses some special challenges. 
Since the TDL frequency can be tuned by adjusting the current and 
temperature, these must be precisely controlled to achieve the 
required frequency stability. The devices are also acoustically 
sensitive and their output is easily amplitude modulated by 
vibrations of any kind. 

To overcome these problems, the TDL was mounted on a 
vibrationally and thermally isolated platform enclosed in an 
evacuated shroud as shown in Figure 7, This platform is located at 
the second stage of a two-stage helium closed-cycle cooler capable 
of operation at ± 10K. A similar mounting scheme is discussed by 
Jennings and Hillman (1977) and has been used aucoessfully to 
stabilize a diode for direct absorption spectroscopic studies. The 
laser and its copper mounting fixture are maintained at cryogenic 
temperatures by coupling the assembly through a flexible braid to 
the second stage of the cooler. This braid, together with a metal 
bellows on the vacuum shroud, isolates the diode mount assembly from 
mechanical shock associated with cooler operation. The entire 
assembly is enclosed in a 50 K radiation shield attached to the 
refrigerator first stage for increased immunity to room temperature 
fluctuations and to eliminate room temperature loading of the 
cooling capacity. A high-precision temperature controller is used to 
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preset the diode mount temperature anywhere between the oooling 
limit of ^9 K, and the maximum recommended diode temperature of 
about 70 K. A temperature feedback loop, consisting of a sensor and 
heater at the diode mount, is used to stabilize the temperature to 
within ±0.3 over a 1 hour period, and has been demonstated to be 
repeatable from day to day. For diode tf8300-l4, this corresponds to 
a frequency stability over 1 hour of about 20 MHz below 1200 cm”^ 
and about 50 MHz between 1200 and 1240 cm”^. Better stability is 
achieved over shorter time periods. 

At constant temperature, the diode output frequency can be 
fine-tuned by varying the injection current, or the frequency may be 
preset to some fixed value for multi-channel integration using an RF 
spectral line receiver. The PbSnTe device now in use exhibits a 
current tuning rate of about 1.0 cm /A near 1180 cm , increasing 
to >\.3.0 om“VA at 1230 cm”\ These rates are alow enough so that 
fluctuations in the frequency due to current supply instabilities 
are negligible. 

Photographs of the refrigerator vacuum shroud, vibration 
isolation assembly and the diode mounting arrangement are shown in 
Figure 8. The top panel (a) shows the vibration isolating bellows 
which externally isolates the cooler second stage from the diode 
mount. Panel (b) shows the diode laser package mounted to the 
copper heat sink. The scale at the left in (b) is in inch units. 

3.3 Beam Matching and Imaging 

The basic requirement for the spectrometer is to combine the 
local oscillator and signal beams using an infrared beamsplitter or 
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Figure 8. Photographs of the diode laser mounting in the cryogenic closed cycle cooler 
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aome other teohnique, and to Image both together at the active area 
of a high apeed photo-mixer which generatea the heterodyne aignal. 

To acoompliah this the ayatem uaea all reflecting optica, with the 
exception of the 50 percent tranamitting ZnSe beam combiner and the 
refrigerator and detector BaFg windowa. The uae of off-axia 
parabolic rairrora (OAP’a) aa a aubatitute for lenaea minimizca 
tranamiaaion loaa and diaperaion effecta and improvea the quality of 
focal imagea throughout the ayatem. In Figure 9, the local 
oacillator path originates with the diode laser output at A, expands 
at f/3 to f/6, ia converted to a parallel beam by the GAP at B, and 
propagates from B to the IR beamsplitter at F, The aignal beam 
follows a similar route, beginning at the instrument focus (I) and 
combining with the local oacillator beam at F. The two beams 
propagate coaxially to the GAP at point G and are imaged together on 
the active area of a high apeed mercury-cadmium- telluride detector 
at H. This portion of the beam combining process alone poses many 
problems, each requiring careful consideration. Gff-axis parabolas 
are difficult to align and produce diffraction limited images only 
in a very small region near the geometric focus. For example, the 
present configuration uses a 2.5 cm beam diameter. The local 
oacillator emitting aperture has dimensions of order 50 microns and 
the detector aperture ia 50-100 microns in diameter. The optimum 
focal lengths in this case are roughly 15 cm and 10 cm for the TDL 
and detector mirrors, respectively. Alignment tolerances for these 
focal lengths are very small and adjustment err ’rs of order 1 mm in 
either the off-axis or axial directions for either mirror will 
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seriously degrade performance. The situation la less critical for 
the long focal length signal beam mirror at J. This is chosen to 
match the f number of the incoming signal beam which varies between 
f/30 and f/120 at the Coude focus of large reflecting telescopes. 

The signal and reference beams must also be parallel and 
coaxial after they are combined at the beamsplitter to avoid serious 
loss of signal. Cohen (1975) has performed calculations of 
heterodyne signal at 10 microns as a function of phase front 
misalignment between the local oscillator and reference beams. This 
study shows that the heterodyne efficiency for matched Airy beams 
will be reduced by 80^ for as little as 4° phasefront misalignment 
between beams. With an f/4 collecting system, this condition is 
reached when one beam is off-axis by no more than several 
millimeters . 

Under normal conditions the diode will simultaneously emit 
several Fabry-Perot modes, separated in frequency by ^ 2 cm”\ In 
order to establish the zero-signal level and perform meaningful 
spectroscopic measurements, only one of these principal modes can be 
allowed to propagate to the detector. To reject the undesired 
modes, an IR grating has been installed in the LO beam at C. The 
grating must allow only one frequency mode to be imaged on the 
active area of the detector. In first order, the reciprocal 
dispersion dv/dp for a grating with element spacing 6 and operating 
angle 3 ia 


d v/d 8 = 6 cos 8 


(3.1) 
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-1 -1 

for V in om . Using tt/2, v^> 1200 cm and a mode separation 

-1 

of <>^2 cm , typical of the TDL used, this leads to mode spacing 
of 


AP 'V 2 X 10"^/ 5 


(3.2) 


produced by the grating. In order to resolve adjacent frequency 
modes, the angular mode spacing Ap must equal or exceed the angular 
width of the grating principle maxima A9 . For a grating illuminated 
by a beam of width D, this is given by 


A9 


\) D cos 3 


. J + « 


(3.3) 


The added term a in Equation 3.3 is the angular width of the diode 
emitting region, and is non-negligible for the dimensions considered 
here. For Ds2.5 cm, fsi5 oA and assuming a 50 micron emitting 
aperture , 


A6 <\i 1.3 3C 10"^ rad (3.4) 

Since AS in Equation 3.2 must equal or exceed A0 the grating element 
density (1/6l must be 
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(1/«S) > 650 cm"'' 


3 65 nun 


(3.5) 


.The requirement guarantees that the modes are resolved but says 
nothing about their separation. In order to have only one mode on 
the detector , the linear mode separation in the focal plane of the 
detector must be larger than the detector diameter. For a detector 
mirror at G with focal length 10 cm, the linear mode spacing AX for 
a 65 line/mm grating is 


aX = f ag 'v 150 pm 


f ^ d \ 


This distance is marginal since it is comparable to the size of the 
detector aperture. A commercially available 150 line/mm grating was 
used which provided a mode spacing of about 350 microns in the 
detector focal plane, isolating single frequency modes on the 
detector. The TDL is mounted so that the polarization of the t 
field vector is horizontal (in the plane of Figure 9) and orthogonal 
to the grating rulings. In this configuration the grating 
efficiency in first order exceeds 85^ at wavelengths between 8 and 
12 microns, and is suitable regardless of the diode used. 

Since the instantaneous heterodyne signal is extremely small 
and masked by the local oscillator noise, a synchronous detection 
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ppooeas is required wbiQd coiaparea the heterodyne aignal piua noiao 
with the noiae in the abaenoe of aignal. In Figure 9f thia ia 
aoQompliahed uaing a rotating blade chopper in combination with a 
atationary flat mirror (0) which providea aignal and reference beama 
from the aouroe in alternate half-oyolea of the chopping cycle. 

Maximum flexibility for a Wide variety of meaaurementa ia 
achieved by meana of dual kinematic mirrora at N. Theae allow a 
calibrated blackbody aouroe L to be placed in either the aignal or 
the reference beam for aenaitivity meaaurementa and inatinanent 
calibration. The blackbody aouroe ia alao uaed for laboratory 
atudiea of aubjeot gaaea by placing a gaa cell in the blackbody 
beam. When obaerving remote aouroea, any remaining teleacope or 
aynchronoua detection offaeta can be removed by changing the aignal 
and reference beam patha uaing a tranalatable diohroic mirror at P. 
The tranamitted viaible beam from thia mirror alao providea a 
conjugate image for teleacope guiding. 

When heterodyning againat the 3un or reference blackbody, the 
aynchronoualy detected aouroe ahot noiae ia appreciable and can be 
comparable in aize to the heterodyne aignal. Thia aouroe of noiae 
ia auppreaaed by ude of an infrared prefilter (Q) to attenuate the 
viaible and near- infrared wavelengtha. Source ahot noise ia created 
since the photocurrent generated in the detector by the source ia 
due to the detector's DC response over all wavelengtha. The 
response of HgCdTe detectors can be quite high, even at near 
Infrared wavelengths where the power per unit frequency is very 
large. An order of magnitude estimate of this effect can be made by 
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comparing the mean square heterodyne signal current from Equation 

$ 

2.10 to the shot noise current from Equation 2,11, with omitted. 

Apart from terms of order unity the relative size of source and 

signal terms is source temperature 'V 5000 K, a 

system atendue of 6,4 X 10**^ and assuming a detector quantum 

efficiency of 'v 0.3 for X > 2.0 pm, the incoherent source power is 

0(30 uW, comparable in size to single mode output powers of the diode 

laser, typically 100-200 yW. Although a large fraction of the 

bolometrio solar flvuc is attenuated by the atmosphere, the 

heterodyne signal for a practical system is typically much smaller 

than that given in Equation 2.10 so that the ratio is still a valid 

indicator of the size of this effect. The IR prefilter installed at 

Q passes wavelengths between 8 and 12 microns. Since the power in 

this wavelength interval is'vll. of the total bolometrio flux, the IR 

filter removes essentially all of the synchronously detected source 

4 

shot noise. Of course, the bolometrio flux varies as T so an IR 
prefilter is not required for observations of sources with 
brightness temperatures significantly less than 1000 K. 

Mention should also be made of the adjustable aperture stop (D) 
located in the local oscillator signal beam. Off-axis beam 
components have been identified as being related to both narrow band 
and broadband noise features observed at IF frequencies while 
heterodyning. The aperture stop is effective in suppressing these, 
as discussed in Chapter V. 

In order to perform remote measurements, the optical system has 
been matched to the 48-inch telescope at the NASA/Goddard optical 
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faclUCy near Greeni)elt, Hd. The incoming beam from the ceiescope 
la approximately f/30 and results in a diffraction limited spot size 
(full width at half maximum) of about 300 mi, crons at the couda focus 
(I) when tho operating wavelength is near 8 microns, The telescope 
beam is collimated using a 30 inch focal length for the signal beam 
parabolic mirror (J) and is combined with the laser local oscillLtor 
beam as previously described. 

3.4 Mathods of IF Signal Processing 

Two methods for processing the IF output of tho HgCdTe detector 
(H) are now being employed. The first method, called sweep 
integration, uses a single channel which is swept in frequency by 
current tuning the TDL. Variations of this approach have been used 
in ptcvioua attemnta to perform diode laser heterodyne spectrcsccpy 
of gas cells and remote sources (Frcrklng and Muehlner, 1977; Ku and 
Spears, 1977; Harward and Hoell, 1980), In Figure 10, the total IF 
signal from the IR detector is amplified 50 dB using a low-noise, 
wideband (0-1.5 GHz) amplifier. A portion of this output is again 
amplified and band limited using a lOQ MHz low pass filter. The 
filter output is square-law detected and the DC signal is delivered 
to a lock-in amplifier and signal averager. Using a high precision 
TDL current supply, the diode laser can be slowly tuned to scan the 
full region of a single TDL mode where diode excess noise is 
minimal. The signal averager is triggered simultaneously with the 
beginn:(ng of each new sweep and permits stacking of both source and 
reference blackbody scans. As a general guideline, in order to 
avoid undue smearing of the line shape, the diode frequency scan 
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IF TO 64 CHANNEL 
FILTER BANK 

Figure 10. Single ehannel, sweep-integration faeteroayne detection electronics 
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rate dv/dt has been governed by the IF filter width fl and the time 
constant t of the look-in amplifier according to 


dv/dt < 0,2 B/t 


( 3 . 7 ) 


For B s 100 MHz, T = 400 msec and a typical single mode tuning range 
of 0.3 cm“\ this corresponds to individual scan times equal to or 
greater than 3 minutes. Since it is scmetimes possible to tune over 
a fraction of a wavenumber, this method' of acquiring data is 
advantageous for observations of highly broadened features, e.g. 
wings of tropospheric lines in solar absorption, or velocity 
broadened features in astrophyaical objects. 

The second approach employs a 64-channel RF spectral line 
receiver (Buhl and Mumma, 1980) to provide multiplexing and extended 
time integration capability. I’his i^s much more useful for 
measurements of doppler broadened features and for sources with low 
brightness temperatures. A block diagram of this system is shown in 
Figure 1 1 . 

In this mode of operation, the TDL local oscillator is pre-set 
to seme fixed frequency in near coincidence with the spectral 
feature of interest, and the IF output of the detector-preamplifier 
is delivered to a variable attenuator at the beginning of the RF 
chain. The IF output is split into 64 contiguous frequency 
channels, each 25 MHz wide, covering 100-1700 MHz. After the signal 
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is divided and passed through an individual filter, it is detected 
and amplified. At this point the receiver noise has been reduced to 
a DC value plus a noise spectrum of several kHz to allow the chopped 
signal to propagate through. Gain pots on the front of the filter 
bank allow the channels to be equalized. The sigpal is then 
integrated in an op-amp circuit for half of a 16 msec chopping 
cycle. At the end of the half-cycle, the individual channels are 
digitized by an A/D converter for each bank and added to the signal 
buffer memory. The integrators are reset (1.2 psec) and the 
integration started for the reference half of the chopping cycle. 
This time the data are stored in the reference half of the memory. 
When the run is finished the quantity (signal - reference) / 
(reference - zero) for each channel is automatically stored on a 
computer disk. The zero is determined by a short integration with 
the RF turned off, yielding individual channel offsets, defined to 
be "zero”. These data can then be manipulated and hard-copied by an 
on-line minicomputer system. Use of a spectral line receiver 
improves the sensitivity of the spectrometer since the entire IF 
spectrum is recorded simultaneously rather than sequentially, and 
long integration times are thus possible. In addition, a 25 MHz 
filter element corresponds to a resolving power u/Av of about 1.5 X 
10^ at 8 microns. This is about a factor of 8 larger than that used 
in the sweep integration mode and very clearly sub-doppler for 
temperatures found in planetary atmospheres. 
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CHAPTSfl IV 

INSTRUMENT DIAGNOSTICS AND TEST RESULTS 

4.1 ■Cnllical.AspeQts of System Performance 

Since lead salt diode laser development is still largely in the 
experimental phase, these devices possess a number of undesirable 
operating characteristics. When used as local oscillators, most 
problems are related to inadequate power output, limited tunaoility 
and a high degree of IF excess noise. Each TDL has its own 
"personality" and must be thoroughly diagnosed to ascertain how it 
will affect the instrument performance. The following sections 
describe the procedures necessary to optimally match the diode 
output to the instrument optics, to suppress the IF noise, to 
determine the diode operating frequencies and to evaluate the 
overall instrument sensitivity. 

4.2 Far-field and Near-field TDL Outpu t Scans 

Optimum matching of the signal and local oscillator fields on 
the detector is essential since only the scalar product . 2,- 
• on the detector active area is effective in producing heterodyne 
signal. Since the TDL output power pattern is often unpredictable 
and Irregular, far-field and near- field scans of the output, shown 
in Figure 12, constitute the first phase of the LO alignment. 

Far-field scans were made as a means of centering the diode 
power distribution on the off-axis parabolic mirror (B) in Figure 9. 
This is necessary to insure that the maximum amount of local 
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(DETECTOR POSITION RELATIVE TO BEAM CENTERLINE - RADIANS) 



Figure 12. Far-field (a) and near- field {"ii) scans of the output 
power pattern from TDL /('6300-lL, operating near 
1185 cm'-^ (T^ = 19 K, = 0.83 ^) . 




oscillator power is delivered to the detector ana that the LO bosm 
is roughly co-axial with the signal beam, as required for optimum 
phasing of the signal and LO fields in the detector focal plane. 

The far- field beam scans were perfonned using a 1 nim pyroelectric 
detector mounted to a 2-axis horizontal/ vertical micrometer drive, 
with the vertical axis motor driven. In Figure 12a the detector was 
positioned about 3 cm from the diode at (A) (less than 2 cm from the 
output window) yielding an angular resolution of about 2°. To 
improve the signal-to-noise ratio, the laser was chopped and the 
detector output and sync signal were fed to u lock-in amplifier. 

The result could then be displayed on a chart recorder. The 
horizontal scan appears roughly gaussian, suggesting a single 
emitting mode, but the multiple lobe distribution in the vertical 
scan is characteristic of high order mode operation. The near field 
pattern was examined by scanning the diode image in two orthogonal 
directions in the HgCdTe detector focal plane (H) , with the grating 
(C) positioned for central order reflection. The image resolution 
is degraded since it is defined by the detector size of >^100 microns 
plus the >/’40 micron blur circle diameter of the detector off-axis 
parabolic mirror (G) . 

Both scans are in qualitative agreement regarding the diode 
output geometry. The horizontal scan shows a gaussian profile for 
both patterns, at least to within the resolution. In the vertical 
direction, however, there is evidence of two emitting sources spaced 
50-100 microns apart, which are just barely resolved. These give 
rise to the fringes observed in the vertical far-field intensity 
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diatribution. The near-field image in Figure 12b was obtained using 
a 10,2 cm focal length mirror, equal tc that at the detector end, 
resulting in a 1;1 reconstruction of the diode image. For longer 
fooal length L.0 mirrors, which seem to yield improved performance 
in the heterodyne mode, some demagnification of the diode image 
occurs and these peaks are completely unresolved. In these oases 
the diode image is positioned by adjusting for a peak in the LO 
induced detector photocurrent. 

The phase relationship between these adjacent peaKs and the 
effects of aperture stopping on heterodyne signal and noise levels 
are not well understood , since these effects vary widely with 
operating conditions and with the device under test. For the TOL 
diagnosed here, the criteria for aperture stopping depend on the 
mode of heterodyne operation considered - the single-channel sweep 
average mode or use of the 64-channel spectral line receiver. It 
appears that off-axis contributions to the LO power are often the 
source of narrow band noise features observable in the IF, similar 
to the spiking resonance symptcm in Figure 4b. Aperture stopping 
provides some discrimination against this out also reduces the 
overall heterodyne signal level. Some IF noise of this type can be 
tolerated using the spectral line receiver since the noise spikes 
generally occupy only one or a few channels and the isolated bad 
channels can be artifically removed from the final data. The single 
channel scheme is much more sensitive to tnis source, since no IF 
discrimination is possible. Attempts made to use RF bandpass 
filters for IF selection have been unsuccessful since the frequencies 
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of these spikes change as the diode is tuned. When this type of 
noise is present, aperture stopping to a diameter of ./'3 cm generally 
gives best results for the spectral line receiver, but a much 
smaller diameter (j'I cm) is gensrally required for a single IF 
channel. This often defeats the sweep integration mode, since only 
a small fraction of the LO power is imaged on the detector. 

^.3 Measurement of Broadband_IF Excess. Nolsfl 

Another dominant source of TDL excess noise is broadband and 
generally covers the entire 0-2 GHz IF spectral range. This is 
generated along with the other characteristic noise types at many 
diode current settings. Figure 13 shows the effect of the mirror 
(B) focal length on broadband noise, examined for two different 
regions of diode operation; one near 1180 cra”^ at a mount 
temperature of >/' 19 K and the other near 1250 cm"' a temperature 
of 42 K. These regions represent the practical extremes of 
operation for this TDL. Curves II and V show quiet diode operation, 
the level being largely due to local oscillator shot noise. Since 
the vertical scale is expressed in dB, the separation i-ii, and also 
iv-v, indicates the ratio of noise power with excess noise relative 
to noise power with the diode blocked. This effect rolls off at 
high intermediate frequencies for ail of the mirrors tried , but is 
minimized for the longer foogi lengths. The LO shot noise 
contribution (ii-iii, and v-vi) is a measure of LO power on the 
detector and hence heterodyne signal, assuming proper matching. 

This also decreases at large focal lengths, simply because mirror 
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(B) no longer captures an appreoiaole portion of the LO power, In 
praotioe the best ail around perfomanoe in the heterodyne mode 
seems to occur using the 6" focal length mirror, In this case, the 
nominal instrument beam diameter collects at about f/4, and captures 
the entire central lobe in the vertical scan of Figure 12a, 
producing the best matching between LO and signal beams, 

4.4 System Sensitivity ■Test U.5ing.J;iie...M-Channel . ■Sp.e.Qtr.al..Llns. 

aeosLy-ar. 

I 

Figure 14 shows a blackbody heterodyne measurement using the 
64-channel spectral line receiver for an LO power of about 300MW, 
and a high quantum efficiency HgCdTe photodiode as mixer. A number 
of different photomixers were used during the instrument development 
and observational phases. This detector was one part of a 12 
element array. Figure 14 shows the overall IF response of mixer, 
preamplifier and cables, typical of LO performance near 1180 cm”^ 
where TDL power is generally high and excess noise minimal. The 
decreased response at large intermediate frequencies results from 
the roll-off of the preamplifer (3 dB B.W, 1500 MHz) and from the 
HgCdTe detector. The standing waves are caused by an impedance 
mismatch between the IR detector and low noise preamplifier. It was 
possible to suppress these by reducing the number of RF fittings 
and the length of cable between the detector and amplifier. 

The vertical scale gives the ratio of heterodyne signal to 
reference which are obtained from the chopped signal and reference 
beams. These values can be compared with the expression for SNR 
from Equation 2.13 (multiplying by 2 for double sideband operation) 
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fco give the system A for these operating oonditions. Values for A 
are shown In tho right hand column. The (snallest &, indicating 
highest efficiency, is about 43 near 500 MHz in tho IF, By analogy 
with the terminology used at millimeter and longer wavelengths, this 
number might be considered to correspond to a system noise 
temperature of 50,000 K. It should be noted, however, that the 
convention for expressing system efficiency in terms of noise 
temperature is not followed for coherent detection near 10 microns, 
For the highest continuum temperatures of interest using this 
technique (i.e. the solar photosphere at </> 5000 K) hv/kT is 0.35, 
so the Rayleigh-Jeans approximation no longer applies, 

4 . 5 Survey of the TPL FabPv-P..e m tJ-lo<le._F_r_eQuensle3. 

The frequency coverage of the spectrometer is limited to 

regions near the diode laser Fabry ~ Perot modes, and the continuous 

tunability is determined by the mode tuning range, In order to 

determine where the powerful modes exist and how widely they tune, a 

* 

survey was performed for diode #8300-14, showing the locations of 
all available modes between 1180 and 1240 om“^ (Figure 15). This 
study is a prerequisite to precise frequency calibration using 
standard calibration lines of reference gases. 

The survey was made by first fixing the diode current at 0,75 A 
and then scanning the mode distribution in frequency with a 
monochromator (S), This was performed at a number of mount 
temperatures over the range of 15K to 45K. While mapping the mode 
distribution, the mode selector grating (C) was placed in central 
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order and acted only as a refleoting element, The LO beam was 
mechanioally chopped using a tuning fork chopper and the sync signal 
was delivered to a look-in amplifier along with that from a 
pyroelectric detector .’.ocated at the monochromator exit slit. This 
is also the arrangement used for making individual frequency 
determinations between measurements during normal operation. 

Absolute frequency calibration of the monochromator to t 0.1 cm""^ 
was achieved by performing a third order polynomial fit to a serie'v 
of HeNe orders, neai’ the first order grating position at 3 microns. 

A few of the sample temperature settings are indicated in 
Figure 15 at the frequencies Mhtr© they produce principal modes. 

The mode locations coincide with the position of the gain envelope 
while their number and size give an indication of the envelope 

shape. The Fabry-Perot frequency separation for this device is 

-1 -1 
roughly 2 cm , with a typical mode tuning range of 0.2-0. 5 cm . 

Examples of continuous tunability can be seen near 1131, 1201, 1207, 

1213 and 1231 where the mode tunes smoothly over the interval 

between two different temperature settings. Tne diode thus lases 

over 15-30^ of its nominal tuning range. As previously discussed, 

however, not even this fraction is usually suitable for heterodyne 

work, since some of it is dominated by large amounts of IF excess 

noise . 

In addition to indicating the mode locations, the survey shows 
the regions in which diode output power is generally high, and also 
the locations where single-mode operation occurs. Although the 
existence of two or more simultaneous Fabry-Perot modes does not 
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imply excesa noiae, the largeat aingle mode power, and hence largea 
heterodyne aignala are obtained during aingle mode operation, 
f!o clear evidence exiata for appreciable output power 
degradation of diode #8300-14 after >/'2 yeara of intermittent uae, at 
leaat in terma of output power, Meaaurementa of long term changea 
in performance are difficult, however, ainoe aeveral different 
HgCdTe detectora have been uaed during the . jvolopment phaae , with a 
wide range of DC and heterodyne reaponaivitiea. Even with the aame 
detector, aome degradation with time oooura, probably at a faater 
rate than for the TDL. 

4,6 Precise Preauenov Callbrat.lon..U.aingLJtef.s.r-£nQ.e-fla5aa 

Although the monochromator is suitable for absolute frequency 

-1 

calibration to within ±0,1 cm , the IF channel width in the 
spectral line receiver is 0,001 cm \ which corresponds roughly to 
the frequency stability of the LO over resonable integration periods 
(>/'1/2 hour). In order to preset the LO for heterodyne measurements 
of selected lines, it is thus essential to determine the LO 
frequency, at leaat to within a few 'milli-wavenumbers . Several 
laboratory gases have been used for this purpose by recording the 
positions of absorption lines that have been identified and 
characterized in published line position studies. The advantages 
and disadvantages of each standard gas are listed below. In any 
case, the initial 0.1 cm”"' uncertainty is small enough to 
unambiguously identify reference gas features . 


%T 
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NgO - line positions itnown extremely well near 1200 
om”^ (Av < 0.001 ora”b but reference lines do 
not always exist over the continuous tuning 
range of a single mode, 

SOg - has an abundance of identified lines, but no 
synoptic observational line posi"..!,;* studies 
for SOg have been made in this jm"?' vf the 
spectrum. The theoretical line positions can 
be in error more than 0.01 cra"\ although the usual 
discrepancy ne£U7 1180 cm“^ is '^.003 cra”\ 
when calibrated against NgO reference lines. 

HgOg - abundance of lines between 1180-1240 cra"^ and 

-1 

positions generally good to ^,003 cm , but the 
liquid is difficult to handle and the vapor reacts 
with BaFg and ZnSe cell windows. 

- generally useful but, like N^O, reference lines are 
infrequent in some regions. Line position 
uncertainties are typically < 0.005 om“\ 

In cases where no N^O line exists near the frequency of 
interest, SOg, and CHj^ have been found to be useful as 

transfer frequency standards. When lines of these molecules can be 
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calibrated against N^O, they provide nearly the same frequency 
positional accuracy as iMgO itself. 

Figure 16 shows an example of a frequency calibration spectrum 

typical of those performed prior to making heterodyne observations. 

-1 

Here, N^O, SO^ and a 0.0625 cm germanium (Ge) calibration etaion 
were used together to achieve a .002 ora uncertainty in LO 
position. The precision was limited largely by the current 
resettability of the TDL output frequency. The calibration spectrum 
was made by alternately placing NgO and SOg gas cells in the LO beam 
while mechanically chopping the LO with a tuning fork chopper and 
plotting the synchronously detected HgCdTe DG output on a cnart 
recorder. The voltage ramp used to sweep the diode current was also 
applied to the chart recorder x-axis drive so that the x-position 
was repeatable with diode current and hence frequency. For this 
measurement the grating was used in first order to select only the 
principal mode of interest. 

The frequency error at any diode current setting is determined 
by the degree of nonlinearity in the diode tuning rate and the 
amount of offset from an N^O line. The quoted SO^ line positions 
are the theoretical values, but the internal consistency for this 
example is quite good (i. .003 cm”^ over the entire page) . The 
fringes from a Ge etaion like the one used here may also show 
appreciable drift due to small changes in room temperature during 
the measurement period. Over the time between a,djacent fringe 
maxima (j' 20 sec), a uniform drift in one direction has been found to 
introduce a 10 to 15 percent change in the absolute frequency scale 


DETECTOR SIGNAL (ARBITARY UNITS) 
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Figure l6. 'SDL frequency calibration spectra obtained using 
H_0 and SO- as frequency standards (frequencies 

‘ -1 ' 
in ca . ) 
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set by the fringe spaoing. 

Also shown on the calibration spectra are the positions of two 
SiO absorption features, the frequencies of which have been 
independently determined by gas phase laboratory measurements (Lovas 
et al., 1981). The spectra of Figure 16 were used to preset the TDL 
current prior to the observations of these SiO lines in the sunspot 
spectrum. The near coincidence of these SiO lines with N^O permits 
an inaependent frequency determination that is in principle as 
accurate as that for NgO, i.e. 0.001 cm”^. In practice, LO 
frequency drift degrades the accuracy of the measurements, 

4.7 Heterodyne Detection of Laboratory M^Q Features 

When using the RF spectral line receiver it is necessary to 
remove the system IF frequency response shown in Figure 14. This is 
done by ratioing the source spectrum to that of a reference 
blackbody. Examples of spectra corrected in this way are shown in 
Figure 17, which shows the heterodyne detection of some laboratory 
NgQ lines against a 1273 K blackbody source. The LO frequency for 
this measurement was determined by locating the N2O feature in 
direct absorption and offsetting the LO a predetermined amount. A 
total of eight 48 second integrations were then made with an N^O 
cell in the reference blackbody beam, followed by an equal 
integration time on the blackbody with the cell removed. Since the 
frequency response of the system is present in both signal and 
reference scans, it can effectively be removed by ratioing the 
signal to the reference, leaving only the gas cell absorption 
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3peo*;rura. The strong and weak fea^:ures are the gv^/RdO) and 3v2~^2 
R(19) lines with rest frequencies of 1177.468 and 1177.550 om . 

For NgO heterodyne measurements the LO frequency can be determined 
to within the line position uncertainty, which is less than 30 MHz, 

_i 

or 0.001 cm at 8 microns. 

In both (a) and (o), the LO frequency lies oetween the two line 
centers placing the observed features in opposite sidebands, which 
then appear folded over each other in this double sideband spectrum. 
The strong feature resides in the lower frequency sideband and the 
weak feature in the upper frequency sideband. In going from (a) to 
(b) the LO frequency has been decreased, the line with the lower 
rest frequency moving closer to the LO line position and the weaker 
line moving in the opposite direction. This figure illustrates the 
importance of continuous LO tunability, not only to accommodate 
large frequency changes, but also to enable determination of the 
sideband in which a given feature resides. The former degree of 
freedom is not readily available with gas laser systems. 

4.8 Comments .Concerning the Remote Detection Tests 

The remote detection capability of the spectrometer has been 
tested in both the sweep average and spectral line receiver modes at 
the NASA/GSFC 48" reflecting telescope, and also at the NASA 3-meter 
IR telesejpe at Mauna Kea, Hawaii. A number of observations were 
made to "shakedown" engineering problems that occur in the remote 
observation mode and verify the field-worthiness of the 
spectrometer. These efforts included; measurements of absorption 
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from terrestrial N^O, CHj^, 0^ ana H^O; oontinuum measurements of 
Mars and Venus (the first such measurements cf Venus); the first 
oontinuum measurements of the infrared object IRO+10216 using this 
technique. These studies are described in more detail in Appendices 
B and C. 


4.9 Summary of Instrument Performance and. Spg.o.lf.lgA.t.Lqjis 
Table 2 lists the specifications for the TDL heterodyne 
spectrometer, as measured using TDL #8300-14. These specifications 
are only representative, and will improve dramatically with the 
substitution of better diode lasers and photodetectors. Some of the 
parameters i\re fixed by the instrument design, while others are 
diode dependent and vary widely, depending on the choice of 
operating frequency. 
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TABLE 2 

Specifications of the Diode Laser Heterodyne Spectrometer 


Instrument! 

Nominal Operating Frequency Range - 
Intermediate Frequency Bandwidth - 
Frequency Resolution ('vJ/AV ) - 
Signal Beam f/number Range - 
TDL Temperature Control Range - 
TDL Injection Current Range - 
Ambient Operating Temperature - 


-1 

800 - 1250 cm 

-1 

± 0.0567 cm 

6 

^ 10 
20 - 120 
9-100 K 

0-2 A 

<32 C 


Diode Laser: 

(Speer, "a-Phys1cs Inc., Laser Analytics 
Div., Diode #8300-14) 


Operating Frequency Range - 

1170 - 1260 

-1 

cm 

Tyolcal Single Mode Tuning Range - 

0.2 - 0.6 

X 

cm 

-1 

Observed Mode Separation - 

2 

cm 

Single Mode Output Power (typical) - 

60 - 300 


Safe Operating Temperature Range - 

< 70 

K 

i 

Overall: 

NEP at 8 Microns: Best - 

-18 

1.1 X 10 

W Hz 

Typical - 

-18 

3 X 10 

W Hz 

Operating Frequency Uncertainty - 

0.001 - 0.010 

-1 
cm ^ 

-1 

cm 

Frequency Drift: < 10 min - 

0.002 

S 1 hr - 

0.005 

••i 

cm 
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CHAPTER V 

HETERODlfNE OBSERVATIONS OF SiO IN SUNSPOTS 

Table 3 lists the dates and relevant data on the series of 

observations of SiO absorption features in sunspot umbrae, made in 

March, 1981. Measurements were made of 5» Vsl-0 fundamental and 

-1 

vs2-1 hot band absorption transitions near 1200 cm (8.3 miorons) 
using the diode laser heterodyne spectrometer described in Chapters 
III and IV. The spectrometer was located at the coude focus of the 
48 inch reflecting telescope at NASA Goddard Space Flight Center, 
Greenbelt, MO. The »5i0 lines were observed at subdoppler resolution 
using the 64-channel spectral line reoeiver. 

The observational conditions were widely variable during 
March, and large differences in air mass, haze and extent of cloud 
cover prevailed during the observations. This variability 
complicates any comparison of observed lines based on their relative 
absorption. Individual umbrae ranged in size from about 5 to 18 
arsec, with values of u between 0.75 and 1.00, As is customary, p = 
cos 0 with 0 the angular distance from solar disk center, With the 
exception of the large spot visible near the center of the disk on 
March 25 and 26, all of the observed sunspot umbrae were 10 arcsec or 
less. The combination of poor observing conditions and small spot 
size made the observations susceptible to photospheric scattered 
light and it was necessary to correct for this effect as described 


! Observational Data 
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in Sec Won 6.2, 

5,2 MiaU&..oLfrMsJi.bmmhiona 

5.2,1 ■Ctoiaa,9t.,QJb,aa£Y.8d. 51^ hims 

Aa shown In Figure 18, the obaarved lines are allowed elecbrio 

-1 

dipole bransitions ab frequencies near 1200 ora wibh a change of ii ‘' 
in bho vibrabional quanbum number v and a change of ±1 in bhe 
robabional quanbum number J, Absorpbion bransibions for which J 
changes by +1 are called R-branch bransibions and bhose for which J 
decreases by 1 are called P-branoh bransibions, A bransibion is 
cypically idenbified using bho nobabion v'v" I (J”) , where bhe 
double primes refer bo bhe lower sbabe and primes refer bo r.he upper 
sbabe , 

The SiO bransibion frequencies used in sebbing bhe local 
osoillabor frequency were bhose reporbed by Lovas eb al. (1981). 
These aubhors have performed gas phase measuremenbs of Si ab 
bhe Nabional Bureau of Sbandards using a bunable diode laser vnd 
have combined bhese resulbs wibh microwave measuremenbs bo yield a 
sob of emperical Dunham ro-vibrabionai oonsbanbs. Their resulbs 
provide bransibion frequencies accurabe bo wibhin t .002 cm for 
bhe 5 SIO lines of Inberesb. The Dunham coefficients obtained by 
Lovas eb al. and the expression used to obtain bhe SiO level 
energies from bhe coefficients are shown ab bhe right in Figure 18. 

The choice of sunspot SiO lines was determined by accidental 
coincidences with the frequencies of the diode laser modes, and also 
by avoiding the regions of heavy waber vapor absorpbion. Figure 19 
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DUNHAM COEFFICIENTS 
(cm"') 


Y|q« 1241.54304 (51) 
Y2o*“S*973907 (263) 
Yjq* 0.0059313 (333) 
Yq|« 21787.4809 (43) 
Y|| "-I5I.0322 (49) 
Ygi* 0.07156 (110) 
Yq2"“29.749(58) 
Y|2"-0.02I9 (38) 


SiO LEVEL ENERGIES 
(cm*' ) 

F-S Y,|(V1--5)* J*( J+l) 

'L; 


LOVAS at al. , 1981 


a 
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Figure IS. SiO energy level diagrua and Dunhoa coefficients 
for determining level and transition energies. 
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isnow3 high resolution ground based FTS measurements by Goldman and 
Blather wick (1980) of terrestrial atmospheric absorption near the 
frequencies of the sunspot SiO lines. These measurements were made 
at low solar zenith angle (about 5 airmasses) to enhance weak 
terrestrial atmospheric features, particularly those of trace 
constituents’. The positions of each of the 5 SiO lines are 
indicated on Figure 19 to show their proximity to major terrestrial 

I , 

absorption features. For some frequencies, tropospheric H^O is 
clearly the dominant observer and observations near the wings of H^Q 
features are very wensitive to the amount of tropospheric water 
vapor. The 1-0 P(20) line in Figure 19c, for example, was sensitive 
to atmospheric haze and one of the P(20) measurements was terminated 
early because of a gradual decline in solar continuum signal during 
the measurement. For the remaining lines, atmospheric absorption is 
not as severe as 'Figure 19 would suggest, since the heterodyne 

observations were made at small airmans (<1.6). 

#» 

5.2.2 Observational Method 

The synchronous detection scheme and observational sequence 
for the SiO measurements was chosen toj remove the terrestrial 
absorption spectrum to reveal only solar features; have separate 
access to both sunspot and photospheric absorption spectra for 
purposes of identifying terrestrial features and; obtain absolute 
sunspot flux estimates. Separate measurements of the umbra and of 
the surrounding photosphere were made by placing a blackened metal 
screen in the reference beam to provide a room temperature reference 


C 
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and by chopping souroo againati aoroen. The aequenoe for each SIO 
obaarvabion oonaldbed of fcha following! 

(i) N Qonaaoubiva paira of loeaaupementia of aunapob againab 
aoraan follow«3<i by phoboaphera againab aoreon. Thia inberlaoing of 
urabral and phoboapherio aoana ininiraiaaa bha affooba of bime 
variabiona of abmoapherio hranamibbanoe and ohangaa in LO 
perforraanoa on bha ooinpoaibo rabio apeobra* 

(ii) M oonaaoubive raeaauretnonba of blaokbody againab card 
followed in aome oaaoa by an N^O cell inaerbed in bha blaokbody 
beam. 

(iii) Poab w meaauremenb oheok of LO frequency drifb by 

checking bha repeababiliby of diode currenb aebbinga ah hhe 

looabiona of reference or S0„ xinea in dirocb abaorpbion. In 

c 2 

bhe cage of bha P(14) maaauremanb, bhe reference NgO line ia close 
enough bo examine in bhe haberodyna mode and sbep (iii) is nob 
necessary, 

Frequanoy oalibrablon prior b6 each raensuremenb was 
accomplished in bwo sbages, using bha mabhod discussed in Sacbion 
4,6, A monoohroraabor was used bo preseb bhe LO Frequency bo wibhin 
'^0,1 om"^ and aibher N.O or SO_ was used as a Vrequanoy sbandard bo 

2 c 

reduce bhe unoorbainby bo several milli-wavemirabers or less. Prior 
bo each raeasuremanb, bha SiO rasb fraquenoias were oorracbed for 
solar robabion ab hhe bime of obsarvabion and posibion on bha solar 
disk. The TDL local osoillabor was bhon adjusbed bo looabo bhe 
feaburo ab bha desired inbermediabe frequency. The 1-0 P(2o) and 
2-1 P(14) SiO bransibions were close enough bo reference NgO line 
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positions to enable heterodyne measurements of both features at the 
same LO setting. Since the NgO 1-in® positions for J" < 60 are known 
to within iO.001 cm"^ in this frequency range, independent 
measurements of these SiO rest frequencies could be made to within 
±0.003 cm"\ the additional uncertainty resulting from small amounts 
of local oscillator frequency drift during the measurement. 

Figure 20 shows a sample data sot for the detection of the 2-1 
P(14) line (t >/' 770 sec), This example is instructive since the SiO 
feature was strong enough to be clearly seen in the un-ratioed 

f 

spectrum (b) and since an absorption feature at v s 1195.695 
cm"^ resides in the opposite sideband. Terrestrial NgO can be seen 
in the solar spectra (a) and (b), as confirmed in (o) where an NgO 
gas cell has been placed in front of a local blackbody source. 




of the photosphe.'’s integrations in (a) and that 


for the sunspot in (b) represent independent measurements, both show 
highly correlated noise features. These result since both 
measurements are ratioed by the same blackbody average, taken about 
1/2 hour later. Over this period of time, small variations in the 
diode laser intermediate frequency noise spectrum occur; enough to 
introduce large amounts of anomalous noise in the ratio spectrum. 

The sequence of collecting data was chosen to enhance the 
detectability of sunspot SiO features at the expense of 
signal- to-noise in the terrestrial absorption speotrura. Because of 
the tight interlacing of sunspot and photosphere scans, and since 
the relative umbral intensity is independent of blackbody, this 
noise virtually disappears in the ratio spectrum (d) which shows the 


ORIGINAt PAQc (S 
OF POOR QUALITY 



Figure 20. Sample sunspot SiO absorption 
measuremant using the TDL 
, heterodyne spectrometer. 
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SIO feature. 

5<3 i^aAMd. of...Fittini!i: Qbaerved Double. Sideband. Features 
The sunspot SiO observations using the TDL heterodyne 
spootrometar are shown in Figure 21. The ordinate for each speotrym 
is the umbral intensity ratioed to that of the surrounding 
photosphere which, to first order, removes the frequency dependence 
of terrestrial absorption. Continuun signal- to- noise ratios for 
these spectra wore typically 300 to 400 near the zero of IF, 
decreasing to less than ^0 at high frequencies because of decreased 
response from the photomixer and IF amplifiers (Figure 14). For 
this reason the absorption features were positioned at IF 
frequencies below HHz. 

I 

From each ^0 observation it is possible to extract 4 
quantities; (i) the continuum intensity and brightness 
temperature, (ii) the line equivalent width, (iii) the line full 
width at half-maximum (FWHM) and (iy) the residual intensity at 
line center. The data in Figure 21, however, show the average of 
the line profiles in one sideband with the continuum in the other. 
Also, in oases where the line center is near the zero of IF, part of 
the line shape extends into the opposite sideband, resulting in a 
distorted double sideband spectrum. In order to extract the 
unfolded line shape fran the data, and estimate quantities (i)-(iv), 
the line fitting routine described in Appendix A was used to 
generate synthetic double sideband gaussian line shapes. These 
profiles were plotted over the data, and good visual agreement was 



Figure 21a. Observed v = 1-0 SiO absorption features 







0.82 
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Figure 21b. Observed v = 2-1 SiO absorption features- 
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obtained by adjusting the line oenter frequency, line depth and 
raioroturb\ilenoe. The doppler width for the SiO iinea is roughly 
constant over a reasonable range of temperatures in the region of 
line formation (Av^/v s 6.0 X 10**^ at 3500 K ana 6.5 X lO"® at 
4000K) , 

5. >4 Kesults of the Measureinenta 

Table 4 lists the measured frequencies, sunspot continua and 
line shapi parameters deduced from the observations in Figure 21. 
The umbral continuum intensities measured relative to the 
photosphere are designated i|i'^ for consistency with the notation in 
Chapter VI. The corresponding line center intensities are 
and the measured equivalent widths are designated W . Limb 
darkening corrections for these measurements are negligible, except 
for the R(4) observation. Even in this case, the correction is 
minor and has been treated in Section 6.3. The measured line 
halfwidths in Table 4 refer to the full width at half maximum. 


SiO Measuranents 
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SiO frequencies measured relative to N^O features using N2O Uns position data by Olson et al. (1981). 
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CHAPTER VI 


ANALYSIS 


^•1 -Thg Llne...Smtheais ProRPam (SPCTRMV 

6.1.1 Qufcllna„Qf the Program 

In order to use the SiO Boasurementa to teat the model 
temperature-praaaure profllea, It ia neoeaaary to have a meana of 
oomputlng line and oontinuua intenaltlea baaed on the Input model, 
Caloulationa of SiO lino and continuum intenaitiea were made by a 
modified veraion of a aynthetio apeotruB program, called SPCTRM, 
4‘iacribed by Deming (1976). Figure 22 ahowa an outline of the 
sequsnee of steps Used by SPCTHM in oraaputing the aunapefc surface 
intenaitiea and lino equivalent widthj. The program accepts aa 
input a variety of atomic and molecular data, along with a set of 
temperature/ pressure points which define the model atmosphere. From 
this it computes the atomic, molecular and electron number 
densities, partition functions and opacities. The radiative 
transfer equation ia then integrated through the atmosphere to yield 
the emergent line and continuum intensities and the line equivalent 
widths. The synthetic line shape information can then be compared 
with the observations, 

6.1.2 Computing the Humber Densities o f each Speoiefl 

The pregram uses subroutines from ATLAS5 (Kuruoz, 1970) to 
obtain the equilibrium nvanber densities for all species (i.e., 
atoms, ions, electrons and molecules) after reading the abundance of 
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Figure 22. Flovchert for the line »ynthe#is program SPCTIU<. 
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all eleoenta relafciv© to hydrogen. Thoae eleraenta are taken to bo 
Independent of depth in the atmoaphere. 

The aiUorn abundance for the SIO oaloulationa waa taken to be 
the photoapherio abundance aa determined by Lwibert and Warner 
(1968). Theae authora obtained a value of Ios^q N(Si) • 7.55 ± 0.10 
on a aoale whore log^QN(H)*i 12.00, uaing obaervationa of aelaoted 
Six linea and Intermediate coupling oaloulationa for the Sil 
apoctrum. The roaulta are oonaiatent with an earlier eatimate of 
log N(Si) s 7.50, derived by Goldberg, Muller and Aller (1980). 

To determine the number denaitiea for all apeoiea, ATLAS5 aeta 
up oonatraint equations involving the numoer conaorvation for each 
element, plus two additional equations; one for charge oonasrvation 
and one for the total density. This results in a 2-diraenaional 
matrix of non-linear equations with a large number of cross terms, 
which are iteratively solved until a self-consistent solution for 
the neutral atomic densities are obtained. The density of the other 
species (i.e,, ions, molecules and electrons) can then be computed 
in straightforward fashion. 

At photospheric temperatures and below, H is largely neutral 
and the main contributors to the ambient electron density are the 
less abundant heavier elements (metals) because of their low 
ionization potentials, The ionization equilibrium for each heavy 
element (i) that donates to the electron density is given by a Saha 
relation of the form 

""l “e ' "i = *1 


( 6 . 1 ) 


whepe 'Jj, (T) ia knovm fijnotion of T involving the ionization energy 

of apeciaa i, All of the heavy elementa thua form a coupled aeb of 

equationa, whoae aimultaneoua aolubion yielda the electron denaity. 

Partition funotiona used in aolving the Saha squationa for moat of 

the atomic apeciea (including H“) were tabulated by Drawin and 

Felenbok (1965) and Moore (1949)* 

The molecular concentrabiona are obtained from the atomic 

number denalties through a aet of temperature dependent equilibrium 

con&tanta and diaaooiation energies* For SiO, it is essential to 

have the ;!«, ■ available estimates for the dissociation energy and 

the equilibrium constant K (T)* The best current estimate for D ia 

n 0 

given by Hulser and Herzberg (1979), with a value of 8.24 ± 0.13 eV. 
Valuea for the SiO equilibriua cenatant K^(7) , defined by 

n(SiO) = K^(T) n(Si) n(0) (6.2) 

were computed from the corresponding pressure equilibrium constants 
Kp(T) over the temperature range 2800K-4000K. These data were 
supplied by S, Abramowitz of the National Bureau of Standards 
(Abramowitz, 1981). Table 5 lists the SiO line strength and 
dissociation data required by ATLAS5. 

A aet of coupled equations of this form also exist for oxygen 
since it is involved in a number of reactions of this type. Those 
included in the model dissociation equilibrium are (in order of 
decreasing Importance); CO, SiO, B^Q, OH, SO, PO, 0^, t-feO, AlO, NaO 


(a) SiO Line Strength Data 
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and CIO. For other molecular apeoies the equillbriuro constants were 
obtained by interpolation ft’otn a compendium of high temperature gas 
thermodynamic properties by McBride at al. 0963), 


Continuous opacity in the umbral atmosphere at the 

wavelengths of interest is due entirely to bound- free and free- free 

absorption by the negative hydrogen ion H". Vernazza at al, (1976) 

have examined the source function for the H" bound- free transition 

in the photosphere at a variety of wavelengths and over a range of 

optical depths which extend into the lower chromosphere. They find 

that wherever H" is an important absorber, no significant difference 

exists between continuum intensities computed using LTC and non-LTF 

bf 

analyses* In LTE the H bound-free opacity is 


, bf, n(H-) 


a^ [1 - exp (-hv/lcT)] 
* 


where 


( 6 . 3 ) 


n(H“) 


.3/2 


2v mUT 




exp (8762/T) 


and a^ is the cross-section taken from Gingerioh (196H). llie H 

f f 

free-free absorption ooeffioient is 
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k/'h n^F/T) (6.4) I 

p j 

i 

j 

The term Equation 6.4 is an analytic fit to the quantum f 

mechanical dipole-length calculations of the H" free- free opacity by 
Stilly and Callaway (1970). 

,1 


F^(T) = 

Included in F^(T) is the correction for stimulated emission, which 
is important at the temperatures and frequencies of interest here. 

The effect of scattering can be shown to be oompletely 
negligible compared to absorption processes. From Vernazza et al. 
(1976), the Thompson scattering absorption coefficient is 

= 6.65 X 10“^^ng (6.5) 

and the absorption coefficient for Rayleigh scattering is 

°Ray ■ "hi 

The cross section R(X) has been computed by Gavrila (196?) for 


^Cl.37 X 10"^^ + (4.375 X 10“^° - 2.60 X 10"'^T)/\)] 
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X < 2.3 miorona and Is a rapidly decreasing function of wavelength. 

The H“ free-free opacity, on the other hand, is rapidly increasing 

with wavelength. As an upper limit to the scattering contribution, 

ff 

the ratio (o., + o ,)/k has been evaluated at X>/'2.3 microns and 

n^y V 

«2 

found bo be < 10 for densities and temperatures typical of the 
umbral atmosphere. The ratio at >/' 8 microns should be several 
orders of magnitude smaller. 

The SiO line opacities at each layer are obtained from the 
transition oscillator strengths lower level number 

densities for SiO (n^„j„) obtained from the total SiO and the usual 
Boltzmann factors. The frequency dependence is given by a 


normalized Voigt profile H(a,v) with v=(v-v^)/AVjj and the pressure 


broadening parameter a. 


Tha exorsaaign for 


e^ 

kJ s ^v'J'v''o«‘-’ " exp(-hu/kT)]ny„j„(SiO) H(a,v) (6.7) 

mcAVp 


The transition oscillator strengths in Equation 6.7 were 

determined using the experimental band transition probabilities 
\tv» listed by Hedelund and Lambert (1972). These values are 
5.11 seo”^ and 9.96 seo"^ for v'-v" = 1-0 and 2-1 respectively 
(Table 5). The transition probabilities are converted to band 
oscillator strengths fy,^,, by the standard relation (Mihalas, 1970; 
Hedelund and Lambert, 1972) 
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^Viv" * (6.d) 

In Equation 6.8, Vy,yi, is the average transition frequency (om”b 
over the band. The band center (Q branch frequency) was used for 
Vyiyii case. All other symbols have the usual meanings. The 

absorption band f values were then converted to absorption 
transition f values t’y,j,y„j„ 8y assigning a lower level 
multiplicity according to whether the tranhition belongs to the R 
branch or the P branch. Following Beer, hambert and Sneden (197^), 

^V'J'v"J” " ^v»v" ^^v' V'^^v' J'v''J“^ 

^ Clm| F(v'v"J'J")/(2J" + 1)] (6.9) 

The ratio |ra|/(2J + 1) is called the Honl-London factor for the 

II II 

transition with m s J +1 for an R-branoh transition and m s -J 

II t 

for a P~branoh transition, 2J + 1 being the total lower level 

. I II I II . 

multiplicity. The term F(v v J J ) is the Hennan-Wallis factor, 

II 

which lies between 0.98 and 1.09 for J,J ^ 78. A value of Fsi has 
been used in this analysis. 

At sunspot temperatures (3500-H000K) the SiO vibrational 
population ratio n( v=i )/n( v=0) is «/’0.6. The larger v=2-1 oscillator 
strengths compensate for the population cfifferenoe , however, so the 
1-0 and 2-1 line intensity ratios are of order unity for equal j. 

The assumption of vibrational equilibrium for SiO was also chocked 
and found bo be valid for the sunspot environment. Computations 
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pepfortaed by Nuth and Donn (1961) for ooiliaiona of SiO with atomio 

hydrogen ahow that for all temparaturea 3000 K tha kinetio and 

vibrational temparaturea ai’o very nearly equal for H preaauirea > 

■»5 

10 atmoapherea. For oorapariaon, SiO formation in aunapota oooura 
at layora whore preaauroa oxoaed 10“^ atmoapherea. 

It is also of interest to examine the terma that make up the 
total level multiplioitiea for tha val-O and va2-1 banda of SiO 
ainoe these direotly determine the lower level populations and hence 
the observed intensities. For a heteronuolear diatomic molecule, 
where interchange symmetry does not apply, the nuclear spin 
contribution to the totnl otatlotlonl welsht of a Invol In 
(2Ig^+1) (2Ijj+1), and being the nuclear spin quantum numbers, 
The nuclear spin part g^^al for SiO since I(^^O) = I(^®SI) a 0 
(Steinfald, 197^). Also, SiO has a simple elaotronio ground 
state configuration. This means that both the electronic orbital 
and spin angular momentum quantum numbers A and E for the electronic 
ground state are zero. As a result' the electronic contribution 
g-=l, and the total level multiplicity, including rotation, 

b 

ti 

S s SpgjjSg s g^ 5 2J +1 for SiO, as expected, 

6.1. A Computing the Surface Intensities and Equivalent Widths 
In differential form the radiative transfer equation at a point 
in the atmosphere is 

p dl^/dz = p I^k^ p (6.10) 


with; 


dz a length element in the vertical direction 
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ii 003 0, Where 0 zenith angle on the sun 
jy the masa emission ooeffioient in erg gro"^ 3eo*'^Hz’^ 
ster“^ 

p the density in gm 

the intensity in erg cm” seo"’ Hz” ster" 

the total absorption ooeffioient, including oontinous 

opacity k^°, line opacity kJ' and scattering o^. 



Since scattering is negligible, only the absorption opacities 
contribute to k^, Defining an element of optical depth dt^ s -k^ p 
dz, and the source function s » equation 6.10 can be 

rewritten 

U dl^ / dT^ s - (6.11) 

The optical depth scale is a direct' measure of the absorption along 
a path which in turn depends on the local densities and temperature, 
for this reason, it i,s a more useful indicator than geometric depth 
in the atmosphere. The source function is in general very 
complicated since it involves the coupled radiation field and the 
bound and continuum states for all particles in the stellar gas. 
Under the assumption of l.i!E, however, absorption and emission are 
related by Kirohoff’ s formula 


j. = ( + ■=/ ) 


( 6 . 12 ) 
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so that 

S^(T) s B^(T) (6,13) 

Using these definitions and assuaptlons, Equation 6.2 oan be 
integrated over optioai depth to got the surface intensity; 


ly (o,u) 


G0 

I" 


(T) exp (-T^ / u) (1 + /x^°) dT^° 


(6.1U) 


with 


T 

V 


V 

! 


(1 + K *• / k °) V 


Finally, the equivalent width for each SiO line profile can be 
evaluated by SPCTRM for direct comparison with the observations. 


W 


V 


J 

Line 


Cl - I^(o,u)/I°(o,u)] dv 


(6.15) 
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In Equii'oion &.15j the integral is taken over the line profile, ana 
I°(o,u) is defined as the surface continuum intensity in the 
vicinity of the line, 

6.2 Method of Correotin g Measured Equivalent Widths 

Before the observed SiO equivalent widths can be used to 
predict the umbral temperatures, they must be oorreoted for the 
effects of s'jattered light, Because of marginal observing 
conditions during most of the SiO observations, the effects of 
scattered light from the surrounding photosphere were found to be 
significant. This was particularly true for the R(4) measurem«;hb , 
with umbral diameters less than «/' 7 arosec, and for the Mar. 17 
P(25) measurement with noticeable telescope wind shake. The results 
suggest an effective heterodyne field of view comparable to or 
greater than the smaller umbral diameters, which is almost entirely 
the result of observing conditions. By comparison, the 8 micron 
diffraction limited field of view (HPBW) for a 48 inch telescope is 
vo 2 arcsec. Contamination from photospherio light is suggested in 
the observational data of Tables 3 and 4, which show a correlation 
between spot size and measured umbral intensity. The small spot 
R(4) observations on Mar. 31, for example, show intensities between 
0.81-0.85. On the other hand, the observations of the 2-1 P(14) and 
1-0 P(26) features on Mar. 25 and 26 and the third 1-0 P(20) 
observation on Mar, 22 were made in larger umbral regions and show 
intensities between 0.73-0.76. This suggests contamination by 
photospherio light in the smaller features. There is also no 
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evidence to indicate that even the largo umbrae fill the heterodyne 
field of view and henoe that oven these lower intensity ratios are 
reliable approximations to the actual umbral continuum intensity. 

It is also not clear whether the relation between spot size and 
measured intensity is totally due to photospherio contamination) 
sinoe changes in continuum brightness temperature of several hundred 
Kelvin between measurements of large and small spots have been 
reported at 2.3 microns, the temperature increasing as the spot size 
decreases (Hall and Noyes, 1969). The 2.3 micron continuum, 
however, is formed where moat models indicate a steep temperature 
g'adient, and it does not follow that the same intrinsic temperature 
fluctuations will be observed at higher layers where the 8 micron 
continuum is formed, and moreover where SiO line formation occurs. 

Figure 23 shows temperature as a function of both optical depth 

and atmospheric pressure P for three competing sunspot temperature 

S 

models, along with previous sunspot continuum observations used as 
boundary conditions for these model's. The supporting observations 
are previous sunspot continuum measurements in the 1-4 micron range 
by Henoux (I 968 ) and Hall (I 970 ), and visible sunspot measurements 
by Maltby and Staveland (1971) corrected for stray light effects 
using tte» Mercury transit at the time of observation. These 
measurements and models are discussed in Section 1.4. All of these 
earlier relative intensity measurements were converted to brightness 
temperature using an empirical photospheric temperature model by 
Vernazza et al. (I 976 ) which shows disk center brightness 
temperatures of 5090 K to 5290 K at 8.3 microns. Supporting 
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CcnparisQn of aunspot continuum observations vlth severaJ, published 
temperature models. The SiO line Center and continuum contribution 
functions are also show. 


1X3 


obs9rvationa for the photosphcrio toaperatures were high altitude 
ground based ffleasurenenta by Kondratyey et al, (1965) and Koutohmy 
and payturaux (1968). 

Tha contribution function for the 1-0 P(26) line center haa 
been included in Figure 23 to indicate the region whore the stronger 
SiO lines are formed and also to define the upper boundary for which 
SlO observationa provide temperature information. The aunapot 
models rapidly diverge at layers shallower than log P * 5.9 » since 
no calibrated long wavelength continuum data are available) and wide 
disagreement exists between them in the region of SIO line 
formation. The SiO tomperature sensitivity through dissociation la 
the key factor in removing the taaparafcura ambiguity in these layers 
orf the atmosphere. 

In Qonbrast) the 8.3 mloron continuum contribution function is 
located close to previous reliable continuum measurements. These 
Indioate what brightness temperatures might be expected at 8 microns 
and allow the data to be corrected 'for photospherio scattered light. 
B7r?rasic agreement occurs between all of the observations and 
models for values of log^Q P^ between 5.3 and 5.5. The region near 
log^Q Pg ''*5.9 represents a particularly good choice for a reference 
layer since the observational uncertainty is small and the 
temperature mod el S) fitted to these and other observationa, are in 
good agreement. A temperature of 389 O K will thus be adopted for 
the layer log^Q P^ s 5,9, Tha divergence of the models at shallower 
layers suggests temperatures between 3550 and 3800 K where the 3.3 
micron continuum is formed. The highest oonoeivable temperature 
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near the peak of the confclnuua contribution function ia about 38 OO 
K, which would occur for a tonperature gradient much smaller than 
that predicted by any of the published models » For comparison the 
lowest continuum intensity observed using the TDL heterodyne system, 
0.74 in Table 3) corresponds to a brightness temperature of 
>/'4000 K, so it la evident that photospherlo scattered light poses a 
problem, 

The scattered light correction procedure may be treated as a 
"filling- factor” problem where the size of the umbra is insufficent 
to fill the heterodyne field of view* The heterodyne field of view 
is defined to have unit area on the solar disk and the spot umbra 
inside the field of view has fractional area A. The unbral filling 
factor A has been determined for each observation by comparing the 
measured umbral intensity to the expected intensity near the region 
of 8,3 micron continuun formation in Figure 23. If l(v), are the 
true umbral line and continuum intensities relative to the 
photosphere, and V(v), are the' measured values contaminated by 
photospher^c light, then the fractional lino depth D(v) * 

1-(4(v)/4'q) and the line equivalent width W can be found from the 
measured values (Appendix C) by 

(B5) 


(E8) 


D (v) 


^ A 


and 


DKv) 


W s 


W 


with the fractional area A given by 


A - (1 - / (1 - 


(E4) 


Adopting T « 3550 - 3800 K as a plausible temperature range for the 
region of 8.3 micron continuum formation, and Including the 200 K 
uncertainty in the photospherlc brightness temperature, this leads 
to values of 0.62 - 0.70 for the umbral intensity i}i^. 

Table 6 lists the values of the filling factor A and the line 
depth and equivalent width for each of the measurements corrected to 
lower intensities. The range of values in Table 6 correspond to the 
range of values of acceptable sunspot continuum intensities. 

6.3 Selection of an Empirical Temperature Model 

The adopted emperical umbral temperature model must satisfy 

two principal criteria. When used to calculate umbral surface 

intensities, it must yield agreement with the previous continuum 

intensity measurements. The second requirement is that it must 

reproduce, within the errors, the measured line equivalent widths, 

for reasonable choices of micro turbulent velocity This establishes 

M 

the temperatures where SiO line formation occurs. 
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Figure 24 shows a set synthetic line profiles computed by 
SPCTRM for different assumed miorotubulent velocities, using three 
published umbral temporaturo models as input (Henoux, 1969; 
Stellmaoher and Wiehr, 1970| Zwaan, 1975). fhe rapid decrease in 
SiO absorption because of dissooia-tion is apparent in going from the 
coolest model (Zwaan, 1975) to the warmest (Henoux, 1969). 

Lino saturation occurs using Zwaans (1975) model atmosphere since 
unity optical depth at line center has moved out co layers where the 
temperature gradient is very small (log P s 3.5 to 4.0). 

To check the value of microturbulent velocity which reproduces 
the data, the observed line widths were averaged and compared with 
the profiles in Figure 22. From Table 4, the measured FWHM values 
in dimensionless units (av/v) averaged 1.01 X 10“^ with a standard 
deviation of about 10^. This average is probably high, since there 
was evidence for frequency drift of the local oscillator during 
several of the measurements. This could amount bo an increase of up 
to 0.3 X lO"^ for individual measurbments that make up this average. 
The value of 1,01 X 10"^ then, is considered an upper limit for the 
observed FWHM values. A more reasonable estimate would probably be 
0.90 X 10"^, which excludes any of the line shapes of Figure 24 
(predicted using Zwaan' s model) and also any large values of 
microturbulent velocity Vp. Reasonable values of appear bo be 
0.5 km/ sec or Itss by comparison with the models of Henoux (I 969 ) 
and Stallmacher and Wiehr (1970). 

Because the spectral lino receiver produces a folded line 
shape, and because of dilution by photospheric light, it was nob 



km/sec 
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possible to "fib" the line profiles oomputed by SPCTflM directly to 
the measurements. The selection of a temperature model was 
accomplished by comparing the measured equivalent widths with those 
computed by SPCTHM for each input model atmosphere. This was most 
conveniently done by computing a curve of growth (COO) for each 
model. The advantages of a model curve of growth for the analysis 
of stellar atmospheres is discussed by Mihalas (1970). For a given 
model atmosphere, the line shape and hence equivalent width of the 
line is completely determined once the transition frequency v , lower 
level excitation energy x and line strength gf are specified. For 
medium and strong features, some estimate of the amount and type of 
broadening must also be mads. At a glance, the curve of growth 
shows the dependence of W on lower state energy, the temperature 
sensitivity of W through the choice! cf models, the character of the 
observed lines - whether weak or strong, and how the uncertainties 
introduced by assianptions regarding broadening and the effects of 
scattered light affect confidence in the choice of an appropriate 
model. 

Figure 25 shows calculated curves of growth for the competing 
sunspot models. The predicted equivalent widths are independent of 
broadening for weak lines (small gf values) but depend on broadening 
for stronger features. Observed equivalent widths corrected using 
the method of &ction 6.2 are also shown and agree nicely with the 
Stellmacher and Wiehr (1970) model for small values of 
microturbulent velocity. This model appears to be the best choice 
for agreement with the SiO observations, Zwaans model yields 
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unacceptably large values of W, which might be expected since it is 
entirely extrapolated for Ios^q Pg X 5*2, On the other hand, 
Henoux's modal predicts equivalent widths which are consistently too 
small. The reason for this discrepancy is nob clear since the 
procedure used by Henoux in his analysis of visible Na D lines 
(Section 1.4) should be fairly temperature sensitive. 

The temperature uncertainty in the region of SiO line formation 
can be inferred from the curves of growth for the Zwaan and Henoux 
models and from the corresponding t<imperature profile.^ in Figure 23. 
These models, which lie above and below the observations in Figure 
25 differ by roughly 0.7 in log (W/v) where observations b-e are 
located. From Figure 23, the same models differ in temperature by >/' 
40Q K throughout the region of SiO line formation. Tno choice of 
acceptable curves of growth which agree with tho data are limited to 
a range of < 0.2 in log (W/v), so to first order, the temperature in 
this region is determined to < (2/7) 400K =.115K. 

Some mention should be made abbut the R(4) observation since 
several aspects of this measurement were different from the others. 
The excitation energy for R(4) is considerably smaller than for the 
other lines, making this measurement particularly interasting. 
However, only a few small sunspots were visible at the date of 
observation, and these were located a larger distance from disk 
center thaf,i the others (m = 0.75). The position of W for the R(4) 
line is somewhat higher on the ordinate of Figure 23 than might be 
expected based on the slope d(log W)/d(log gf) of the Stellmacher 
and Wiehr (1970) curve of growth. On© might expect the R(4) result 
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to be located below the curve, ainoe any increase in temperature in 
the smaller spots should increase the SiO dissociation, decreasing 
the equivalent width. One possible source of this discrepancy is 
the u dependence, which has been taken into account in Figure 23. 
This affects the equivalent width directly and also indirectly since 
limb darKening alters the measured continuum intensities, affecting 
the scattered light correction. These effects were checked by 
computing the Stellmacher and Wiehr (1970) and HSRA (1971) models 
over the ran'ge from u=0.6 to 1.0. It was found that relative. to 
disk center both <t>Q and log W inoress^ed by ./’0.02, These numbers 
were used to normalize the observations to disk center, resulting in 
a reduction of This correction is in the right direction, but 

the measured range of W for R(4) remains anomalously high. 
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VII 

SUt^MARI AND CONCLUSIONS 

7*1 ^iaiDa£y-^nd-JLQnglU3lon5..from the Sunspot SIQ 

Apalxaia 

The measurements discussed in Chapters V and VI constitute the 
first reported detection of SiO in the solar atmosphere and the 
first such measurement using coherent detection. Although there are 
observational uncertainties in the measured equivalent widths, the 
absence of blending from adjacent lines in the IR and the high 
temperature sensitivity of the SiO abundance through dissociation 
permit an unambiguous temperature determination over an extended 
region in the umbral atmosphere, By combining the TDL heterodyne 
SiO and continuum observations with previous IR and visible 
continuum measurements, the vertical temperature structure in the 
umbral atmosphere has been aounded with an accuracy of ^ 100 K over 
the region log P = 3.5-5. 5. This permits an unambiguous selection 

O * 

of a temperature model, provided it is well behaved. The published 
umbral model of Stellmacner and Wiehr (1970) yields good agreement 
v;ith the observed SiO equivalent widths, so this model has been 
adopted as the best choice consistent with the TDL heterodyne 
results. From the temperature sensitivity of the SiO abundance, the 
lack of large amounts of scatter in equivalent widths observed from 
spots of different sizes suggests a remarkable consistency of 
temperature in the region of SiO line formation (AT < 80 K) . 

From a comparison of the measured and computed line widths it 
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does not appear neoessary to invoke any broadening processes otner 

than microturbulenoe to explain the Observations, In fact, tne 

lowest reasonable estimate suggested by Steilmaoher and Wiehr (1970) 

of 0,5 km/seo yields line widths comparable to the observations 

when their model is used as input to the line synthesis routine 

SPCTRM, Comparing the TDL heterodyne results with figure 24 

suggests that may be even smaller than 0.5 kra/seo in the region 

of SiO line formation (log P s 3,5 to 4,0). The absence of 

s 

evidence for pressure broadening of SiO in the ambral atmosphere is 
of interest since no laboratory measurements have been made of SiO 
broadening by foreign gases. An upper limit for broadening of SiO 
by the stellar gas oan be obtained by setting V^sO and adjusting the 
damping constant in the model until agreement with the observations 
is obtained. This was done for the 2-1 P(14) line using the 

-4 

Stellmacher and Wiehr (1970) model and yielded a value of 1 X 10 
om"VTorr as an upper limit for the broadening coefficient for SiO 
in the umbral atmosphere, ' 

Regarding future sunspot observations, a number of di- and 
tri-atomic molecules expected in fair abundance in sunspots have not 
yet been observed and some of these have rotation-vibration bands in 
the 8-12 micron region, accessible to the TDL heterodyne technique. 
From Table 1, those species are; VO, TO, PO, AlO, CO^, SO, HgS, Hg, 
NaH, NH^, NS, HON and CS. These species were predicted by Gaur et 
al.(1973) to be formed at both larger and smaller depths than SiO so 
that similar observations and analyses of these species should yield 
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a *:horough characterization of tho umbral uppor atinospharo. 

7.2 Additional 'rPI. Hatgrodyne Qbaervatlonal Problem8_of 

.Qsnoral-Intgr.ggt, 

The sensitivity of the diode laser heterodyne technique ana its 
practicality as a multi-disciplinary observational tool continue to 
improve. Ihese advances are largely the result of recent 
improvements in diode laser fabrication and an increased 
understanding by the user of the noise properties of these devices, 
With each improvement in heterodyne sensitivity, new observational 
problems can be considered which lend themselves to the high spatial 
and spectral resolution as well as the inherent tunability. 

Besides the obvious advantages of diode laser heterodyne 
apectrescopy for the study of molecular lines in the sunspot 
environment and for sub-doppler observations of telluric atmospheric 
constituents in solar absorption, the method has great potential for 
further astrophysical work, Several interesting problem areas 
should be mentioned here; 

1) Anomalous emission features in the solar spectrum, near 
12,33 and 12,22 p have been discovered by Muroray et al, (1981), and 
have been conjectured to originate in the solar chromosphere or 
corona. At the observational resolution of 0.02 cm”\ the features 
were unresolved and exhibited a peak intensity >/’10jS above the 
continuum. The fully resolved line shapes obtainable using the 
heterodyne technique will permit accurate measurements of the 
kinetic and brightness temperatures for these features. In 
addition, accurate determination of the line center positions should 
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provide informaSion on the dynaraica in the regiona where the 
featurea are foraed. 

2 \a demonatrated by She obaorvaSiona of itora, Venua and 
IflC+10216, TDL heterodyne apeoSroacopy ia evolving to the point 
where remote aourooa other than the Sun can be detected with 
aignal-to-noiije adequate for apectroacopic purpoaaa. The meaaured A 
of 43, which occurred uaing the HF apsotral line receiver ia within 
a factor of 7 of the beat performance from the flASA COg heterocyne 
ayatsm. Part of thia difference (a factor of v<2) ia accounted for 
by the reflection charactaristica of the ZnSe beamaplitter , which 
differa between the two ayatema, 

3) Several bright IR atara, e.g. IRC + 10216 (already 
observed using the TDL heterodyne apeotrometer) , VT Cania Majoria, W 
Hydra, and a Orionia have also been considered for study after some 
additional improvement in system sensitivity. A variety of 
molecular absorption and emission lines have been observed in these 
and similar objects at both IR and millimeter wavelengths, The list 
of such detections has grown too large to discuss in this context. 

Of great interest is observed SiO maser emission in a number of 
these IR sources. Infrared measurements of vibrational transitions 
between the levels of interest may answer some long standing 
questions regarding the mechanisms for pumping these masers; 

4) Molecular hydrogen emission has been detected in the 
Orion region (Beckwith et al., 1978 | Beck et al., 1979) and 
observations of selected Hg lines at high spatial and spectral 
resolution should reveal new information about this complex region. 
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Suoh neaaure.T}Qnts8 will be possible for tsne Hg line as 12, 2d aicrons 
when She sysSsra senaltivity is improved and will beoorae possible for 
Hg lines at 17.03 and 28,22 miorona whan long wavalengSh infrared 
mixers and TDL*s become available, 

5) Infrared measurements of some oorapaot HII regions reveal 
large fluxes in She fine structure lines CArlll], CSIVJ , and CHell] 
at 8,99» 10,51 and 12.81 microns (Willner, 1977) # for a 3 meter 
?r?rtelescope aperture, estimates yeild brightness temperatures in 
these lines as large as 300K, making them attractive sources for 
diode laser work, Low resolution observations of these lines have 
also been made in planetary nebula (Aitken et al,, 19T9) ana suggest 
line brightnesses adequate for high resolution observations with 
signal to noise ratios larger than 50. 

7,3 F.dt,ucs,Di£g,sii,2Qa.in..IgL H. 9 t.Si:b!i::n,aJ;tiaA.nfflien,t,, D,fiv.aiaBinanii 

In order to make diode laser heterodyne detection fully 
competitive with (or superior to) ^the overall performances of gas 
laser instruments, advances are necessary in the areas of TDL and 
photoraixer development as well as in the overall instrument design. 

In order to successfully perform most of the observations 
above, it is essential to have TDL's with improved (i) power 
output, (ii) mode purity, (iii) excess noise suppression, (iv) 
tunability and (v) operating temperature range. One additional 
long-standing proolem - that of device deterioration by slow contac 
degradation - has been largely solved. This results in 
substantially increased lifetime over devices of just a few years 
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ago, Itesna (ii) mi (iU) are eloaely rolafesd ainoe Many of *^no 
diode noioe proportiea are coupled to beat effeota and awitoning 
between competing frequency nodea, The elimination of multi-moding 
in atripe geometry laaera alao greatly enhances aingle mode 
tunability by reatrioting operation to only the lowest order 
tranavarae mode. Aa a result, auooeaa in Improving (ii) ehoulc also 
largely eliminate problems associated with (lii) and (iv), 

The HgCdTe IR photoraixors incorporated into the system have 
been provided by D. Speara nf MIT Lincoln Laba (Speara, 1977)* 

These devices exhibit high quantium efficiencies in ahot-noise 
limited applications, out have bean optimized for use with COg gaa 
laser heterodyne aystema which routinely deliver ImW to the 
detector. Because of the inherently low power outputs of present 
TDL's, detectors are required which are optimized for smaller local 
oscillator powers. Such devices are juat now beocming available 
(Speara, I 98 I) and should improve the heterodyne efficiency in cases 
where the available TDL power la s^plit among several usable modes. 
One major engineering problem associated with the spectrometer 
has been the inability to stabilize the diode laser frequency over 
reasonably long integration times, A related problem is the 
frequency uncertainty that arises when it's necessary to extrapolate 
the diode frequency away from precise reference gas line positions. 
In practice, the frequency stability of the TDL should be less than 
the single IF channel width of Z5 MHz over a period of 1-2 hours, 
which is the maximum practical integration time for most sources. 
Present commercially available temperature controllers employing a 
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neater-sensop servo loop are not capable of conaiatently satisfySnj 
this requirement, particularly when the riiodo tuning rate exceeds 

«r1 

</> 2,0 cm 7K, Also, because of non-linearities in the diode tuning 
rate ana because of thermal drift of the Ge etalon fringes, local 
oscillator frequencies normally cannot be accurately extrapolated 

more than >/'0.01 cra“^ from a reference gas line, A reliable method 

*1 

of predicting and holding the TDL output frequency to <0,001 era is 
necessary for reliable oosarvational work and long integration 

times, Poultney et al. (I9d0) have achieved TDL frequency control 

-4-1 - 1 

to v'lO cm , using reference gas features and a 0.01 cm air gap 

etalon. The TDL frequency was determined relative to the etalon 

fringes using an interpolation scheme based on a TDL fm analysis of 

the etalon fringe pattern. Another approach (Savage and Augeri , 

i960) uses a servo-controlled Fabry-Perot etalon in which the 

transmitted LO power is synchronously detected and used to fine tune 

the TDL current controller , A number of other schemes employing 

reference gas lines and/or variations of these methods might be 

employed to achieve the desired fi'*equency precision and long-term 

stability. 

Another design problem results from the inefficiency introduced 
by the zinc selenide beamsplitter. In the present spectrometer 
design, half of the useful signal and LO oscillator power is lost at 
the beamsplitter because of its 50/50 reflection/ transmission 
characteristics. Two plausible schemes for recovering this lost 
signal are: ( i) introduction of a second HgCdTe detector, the 

output of which is combined in the proper phase with the first 
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detector usitvj a broad-oand If coupler » (il) fiepiaoement of the 
beamaplittar with an all-rfifieative beam combiner. This second 
scheme has the potential of being easy to implement, and its 
inclusion, into the spectrometer optical train is presently being 
considered , 
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APPENDIX A: 

TDL riETERODXNE DETECTION OF ATMOSPHEKIC NgO, CH^, 0^ AND HgO 

A number of obeervations of terrestrial features have Dean made 

in solar ocoulation while besting the instrument at the NA3A/GSFC 48" 

reflecting telescope, near Greenbeit, MD. Figures 26 and 27 show 

heterodyne detection of atmospheric absorption features measured in 

the sweep integration mode. Figure 26 shows P and R branch lines in 

the 2vg fundamental band of NgO, The features in Figure 27 belong to 

the V. fundamental band of ozone, based on line identification 

studies and position measurements by Goldman and BlatherwioK (1980) 

and Barbe et al. (1981). Each solar measurement was followed by a 

reference line measurement using an NgO cell in front of a 1273 K 

blackbody. A calibrated etalon and a calibrated blackbody continuum 

spectrum (not shown) were measured for each line allowing 

quantitative analysis of the line shape. Because toe strong NgO 

* - 1 -1 

reference lines have separation larger than 0.5 cm , the 0.1 cm 
frequency calibration using the Infrared monochromator maKes it 
possible to identify each feature and obtain its frequency using NgO 
line position standards. A recent compilation of NgU absorption line 
standards, accurate to < 10 HHz (Olson et al., 1981) has been used 
for this purpose. 

Both NgO scans in Figure 26 were obtained in early June 1980 at 
nearly the same air mass. However, a substantial amount of water 
vapor was present on the day the R(19) line was obtained and 
absorption in the wing of a highly broadened HgO line centered at 
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Figure 27. Heterodyne neasurecent of atmospheric 
ozone lines In near coincidence with 
laboratory integration 
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tl87.0 om’^ accounts for *;he reauced solar continuum signal during 

that measurement, The R(19) result was obtained by stacking 3 scans, 

each with a time constant of 1,25 sec for a total integration time 

(t) of 3.75 seconds. Only a single scan with t = 1,25 sec was 

required for the P(12) measurement because of nhe larger solar 

continuum signal and lower TDL noise level. 

Most of the noise in figure 26 results from the choice of 

amplitude resolution used when storing the scan average in the signal 

averager. The excess noise, however, results from mode competition 

in the diode laser which produces large audio frequency fluctations. 

In practice it is sometimes possible to adjust the diode current and 

mount temperature to move the interesting feature into a region where 

noise is not observed. For the P(12) line the double side-band SNR 

for the line wing is as high as 200 for an IF bandwidth of 100 MHz, 

integration time of 1.25 sec and blaokbody temperature of 1273 K. 

The expected SNR for an ideal system under these same conditions is 

* 

about 16,000, making the measured system delta (A) about 80 in this 
region of the scan, A delta as low as 60 has been observed in the 
single channel sweep mode, although not in close proximity to an 
absorption feature in any of the subject gases used here. 

Figure 28 shows atmospheric methane and ozone observed in solar 
absorption after interfacing the spectrometer to the 64-channei 
spectral line receiver, The detection of CH^ in (a) was confirmed by 
heterodyning a CHj^ cell in front of the reference blackbody 
immediately after the atmospheric measurement. Line position 
standards for methane (Reisfeid et al., 1973) were used to identify 
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the feature, after frequency calibration to within iO.1 cm“^ using a 
fflonochromator , The poor signal to noise in this taeasuretnent results 
from' both excess noise in the diode mode and low atmospheric 
transmission ('^25iO at this frequency, as given by the ratio of 
observed solar signal to reference blacKbocy signal. The CH^ feature 
lies well within the wing of a highly absorbing HgO line at 1198,2 
cm“\ This would account for the heavy absorption, since the 
atmospheric temperature and relative humidity were both high at the 
time of measurement, 

-1 

Absorption from H^O at 117*1.5 cm also limits the signal 
observed in (b). This shows atmospheric 0^ line identifications. that 
were made using 30^ as a frequency standard. The LO frequencies for 
these mea.surements wore determined by direct absorption using a 10 cm 
cell filled with 5 torr of SO^. By noting the LO current settings at 
several SOg line centers and using line position data for the band 
of SOg (Allario et al., 1975) the LO frequency could be found by 
interpolation. Based on SOg abundance estimates (Aikin, 1979), no 
observable SO^ absorption is expected in the solar transmittance 
spectrum. Both of the features in (b) were identified as belonging 
to the fundamental band of 0^ using published ozone line positions 
by Barbe et al. (1981) . 
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APPENDIX B! 

HETERODXNE HEASUREMENTS OF MARS AND VENUS 

Obaervations of Mara and Venus uaing the TDL heterodyne 
apectrometer were made during November i 19^0 at the ooude focua of 
the NASA 3-meter IR telescope at Mauna Kea, Hawaii, A photo of the 
spectrometer in the 3-meter coude room is shown in Figure 29. The 
beam of this telescope is approximately f/120, yielding a 0.5 
arcsec/mm field at the coude focus. Th<' d micron "seeing" at Mauna 
Kea is about 1-2 arcsec under good conditions, and because of the 
altitude (>^4.3 km) the amount of precipital water vapor is very small 
(typically 1-2 mm). Since HgO is the dominant absorber at the 
frequencies of interest, this altitude is of great advantage for IH 
observations. 

The planetary measurements (Figure 30) were made using the 
64”Channel spectral line receiver and the EO was adjusted to operate 
using a single powerful Fabry-Perol; mode. The mode frequency was 
measured with the monochromator to be 1181.4 cm"*^. Examination of 
0.06 cm“^ resolution ground-based atmospheric transmittance data by 

Goldman and Blatherwick (1980) shows that this frequency lies between 

_1 

terrestrial NgO features at 1180.9 and 1181.8 cm . No other major 
terrestrial’ absorption features exist near the observing frequency so 
nothing was expected to appear in the spectra, given the observed 
signal to noise ratio. The reference 1300 K blackbody measuremenr, 
was made prior to the planetary measurements and snows an operating 
system A of 45 near 400 MHz IF. The blackbody measurement is shown 
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in Section 4.4, Figure 14. The synchronous detection aoneffie involved 
the use of a wobbling secondary with a j'I arctain throw, in place of 
the rotating blade chopper. Either beam could be centered on the 
planet with the other on the sky. This oapabllity allows the 
cancellation of systematic differences between the two beams that are 
comparable in size to the planetary signal. 

The Mars result (top) consisted of 10 co-addcd integrations at 
disk center made by alternating beams, for a total integration time 
of i/'64 minutes. The .-lars measurement was made during daylight hours 
since Mars transit occurred at about 2 p.ra. local time. This made 
tracking a major problem since the high background light level made 
the disk of Mars almost invisible in the visual field. Tracking 
difficulties were not quite as severe for the Venus measurement, 
because of the planets higher visible albedo and the high solar 
zenith angle at the time of measurement. 

Adopting a disk center continuum temperature of 260 K for Mars 
and 240 K for Venus, the computed fluxes of these sources when 
ratioed to the blacKbody should be 3.1 x lo"^ and 1.8 x lO”^, 
respectively, The Venus continuum signal of 0.64 x 10“^ then, 
suggests a total transmittance for atmosphere and telescope of >/’0.4. 
Compared to Venus, the ratio Mars/BB should have been «/'1.2 x 10“^, 

The observed value of 0.8 x 10 is presumably the result of guiding 
difficulties. The integration time for Mars was only about naif that 
for Venus, which accounts for the poorer signal to noise observed in 
the Mars measurement . 
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HETERODlfNE HEASUaEMEtJTS OF THE INFRARED OBJECT IRC+10216 

Meaauremonta of IRC4-10216 were alao aiade afc i*iauna Kea during 
January 19B1, Aa Figure 31 showa, thoae meaaureajonba were made on 
January 14 and 15 uaing diode laser modea ab 1219*30 and 1227*12 
am"\ These frequenoioa were eafcaoliahed by first ooarae oalibrabing 
using bhe nionochrcmaber , and bnan of fsebbin? from Ng® referenos lines 
ab 1219.260 and 1227.099 cm"’*. A 1300K reference blacKbcdy was used 
for these measuremenba, 

Aa with the planetary run, aucoeaaive measurements on the source 
were made by alternating beams from the vjobbling secondary to remove 
thermal background offsets that exist between the two beams. In 
order to successfully track and remain centered on the source, the 
brightest part of which is confined bo a region less than ^<1 arosoo , 
the instrument was periodically used as a broadband detector by 
spiralling the telescope over a l-3Varosec field, while sensing for a 
peak in the synchronously detected direct signal from the source. 

This procedure was repeated after each pair of source integrations 
once every 10 minutes) and was necessary to correct for slow drifts 
in telescope pointing. 

One motivation for making these measurements was the possibility 
of detecting SiO absorption in this source. For this reason the 
choice of LO frequencies was restricted to close coincidence with 
expected positions of low J, v = 1-0 and v s 2-1 lines in the source. 
The mode used on January 14 exhibited excessive amounts of excess 
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Figure 31. Heterodyne measyrements of the Infrared 
object IRC+10216, performed at Mauna Kea. 
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ncise an high IF, causing a sheep roil-cfT in hhe IR response. 

Ideally any IF response charaoheristics will disappear in the ratio 
spectra. During the time intervals between source, moon and BB 
reference measurements, however, small changes in LO performance 
occurred. Because of the low signal levels being detected, these 
changes introduced large anomalies in the ratio spectra. For this 
diode laser mode the general IF noise level showed time variability 
as exhibited in both moon calibration runs and in several sky versus 
sky comparisons. 

The January 15 data appear considerably better, partly because 
the TDL mode in this case was well behaved, and partly cecause of 
improved seeing - an important requirement for observations of point 
sources. The heterodyne calibration measurements were made by 
placin' ?ne beam of the wobbling secondary on the lunar limb. Since 
the lunar phase was near first quarter, the subsolar point was also 
near the limb and the lunar continuum temperature was expected to be 
about 39 OK. This would yield a mco^i/BB ratio of 0,032. The measured 
ratio of 0,012 - 0.015 suggest an atmosphere - telescope 
transmission factor of about 0.4, in agreement with the November, 

I 9 S 0 results. The ratio IRC+102l6/racon was measured to be v 0.20. 
This result is not inconsistent with a value for IriC+102l6/mccn of 
0.38 measured by Kostiuk and cowcrkers (1981) with the Goddard 10 
micron COg heterodyne spectrometer at Kitt Peak (Kostiuk, 1981). The 
moon was full at the time of these measurements (January 2, 1980) and 
because of the lower temperature at the lunar limb, larger values of 
the star/raocn flux ratio would be expected. The diode laser 
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measureoienl; suggests an 8 mior?n brightneas hemperanure cf 290K fcr 
IRCJ+10216. 




APPENDIX D: 

DETAILS OF THE VOIGT PROFILE DOUBLE-SIDEBAND FITTING ROUTINE 

The spectral line receiver produces a folcec double sideband 

spectrum which complicates any analysis of an absorption feature in 

one of the sidebands. For two special cases the unfoldins process in 

simplified: If an isolated narrow feature is observed centered at 

some distance out in the IF, e.g., a doppler-broadened line at SAVj^ 

from the local oscillator frequency, then the effect of folding is to 

double the continuum intensity relative to the line depth and reduce 

the contrast in the. line by a factor of 2. The single sideband 

iineshape tnen be retrieved by doubling to get 

2(<> ->» )/^ . The other special case arises when the local oscillator 

is placed exactly at the line center position. Then, for a symmetric 

line, the folding of each half of the line around the LO position 

doubles the line contrast so the measured (41 -4' )/‘t is the single 

c u c 

sideoand value. ' 

In ail other cases, either because of excessive broadening or 
close proximity of the line center to the LO position, folding of the 
line wing in the opposite sideband causes a distortion of the line 
shape and a visibility intermediate to that of the two special cases 
above. This is often the case v;itn the observed SiO features. 

In order to recover the line shape and extract reliable values 
for the widths, equivalent widths and residual line intensities, a 
computer program was v/ritten which folds a Voigt function about some 
selected point and overlays tne folded curve on the observed 
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double-sideband fea*vure. Figure 32 and bhe following discussion 

illustrate the details of the folding and fitting process. 

Letting 4 > and be the continuum and line center 

c 0 

S S D 

intensities, the single sideband intensity 'J>(Av) at a point in an 
absorption line a distance Av from line center is 

s « - [« - )] [H(a,v)/H(a,0)] (D1) 

c c o 

In Equation D1, H(a,v) is the Voigt function, with the customary 
definitions for a and v. 

a = = ratio of ccilisionai broadening 

to doppler broadening 

Here Av^^ is the doppler half-width for the absorption line. Since a 
is dimensionless, y has the units of Av^^ and represents the 
ccilisionai broadening contributiorf. Also, letting v be the line 

w * 

center frequency 

V = av/av_, = (v - v^)/Av,- (L2) 

As usual, V is the independent variable which gives the distance from 

line center Av, in units of the doppler width. Of course v = 0 when 

V = V . Placing H(a,0) in the denominator of Equation D1 normalizes 
c 

SSl3 

the height of the Voigt function at line center to (^i - ^(v )) sc 

C 0 

®^*^4i(Av=o) s as required. The double sideband 


that 
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infconsi^.y at an intermediate frequency is tne average cf the 
single-sideband values at and - v^^) , Frcm Figure 

32 this can be written 


♦'V = 


1/2 


IF 


clF 


ssb 


<t>( V 


IF 


'’clF^ 


(u3) 


where 


ssb 


«(v 


IF 


V = 
cIF^ 


ssb 


♦'W’ 




In Equation D3 we have defined 

H'(a,v^) = H(a,v^.)/H(a,0) 


with 


= (V 


IF 


- '’oIF^^^'^D 


Simplifying, Equation D3 becomes 
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♦(V 


♦. - 1/2 (*, 


SSD 




CH'(a,v ) 


+ H' (a.v ) 


(D4) 


It is also important to mention the treatment of the doppler 

term in this analysis, For a moleoule of mass M (in amu) at 

temperature T, the thermal doppler (1/e) width for absorption (or 

emission) at frequenoy v , hereafter called , is 

c r* n 

= (2 k T/M o^) (D5) 

f 

with k = Boltzmann's constant and M ra^, the molecular weignt (M m„s44 

n n 

for SiO) , For sunspot temperatures (T s 3500-4000 K) Av. for SiO is 

tn 

rather v/ell defined, and lies between .0038 ana .0041 cra“^. The 
effect of microturbulence broadening is taken into account by adding 
an assumed microturbulent velocity in quadrature to the thermal 
doppler velocity to obtain an effective doppler width The total 

particle velocity is related to Av^^ by 


and can be expressed as the sum of the thermal part (5,.) and the 

th 

microturbulence contribution (5^^) by 


r 2 2. 1/2 

- ^^th ) 


(D7) 
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2 2 

;Since 5^.^^ s ■ ( c/v^) Av , Equation D7 becomes 

AV„ = v„/o = (CS) 

From Equation d 6, the themal term in brackets is only a function of 

temperature, Since is usually given in wn/sec, the first terra can 

2 

be expressed in (km/ sec) so that for SiO, Equation 08 becomes 

= (lO^v^/c) [3.76 X 10“^T + (D9) 

The procedure used in fitting the observed lines was to adjust 
sst) 

♦ a and until good agreement was obtained between 

the cbsKi’vaticns and the folded Voigt lineshape (Equation 04), whicn 
is overlaid with the data. The line width and equivalent width could 
then be extracted from the unfolded line shape. 
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APPENDIX E: 


EFFECTS OF SCATTERED PHOTOSPHERIC LIGHT ON MEASURED LINE SHAPES 

For observations of small spots or in cases vjhere the heterodyne 
field of view does not reside in the spot umbra for the entire 
integration period, some dilution of the umbral spectrum occurs 
because of contamination by surrounding photospheric light. This 
situation is equivalent to a filling factor problem where the source 
size is smaller than the field of view. For a beam area of unit size 
with a fraction A of the beam filled by the spot, the measured umbral 
intensity I'(v) will be 


I'(v) = A I(v) + (1-A) 


(El) 


Here I(v) is the corresponding intensity for A=1, (i.e. When the 
umbra completely fills the beam). The measured quantity (^'v) is the 
ratio of spot to adjacent photosphere intensity I'Cv)/!^^ so that 
most useful form of Equation El is obtained by dividing each term by 
I , to get 


i)i'(v) - A 4>( v) + (1-A) (E2) 

Likewise for the continuum measurement 

= (1-A) 

c 0 


(S3) 
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F'iP each measurement, the effective filling factor A can be deduced 
from the measured continuum i)»' after a suitable value for the "true" 

Q 

relative umbral intensity has been determined, Kearranging 
Equation E3 gives 


A = (l-ij* )/(!-* ) (E'^) 

c 0 

For each measurement it is then possible to get an upper limit D(v) 
to the line depth D'(v) at any point on the line, Since 

D*(v) s (it>* “ s A/i)i* (i(i 

0 c c c 

solving for D(v) gives 

D(v) = (* - <^(v))/t„ s (^' A) D<(v) (E5) 

C C w w 


The equivalent width W can be adjusted in a similar way. From the 
definition of equivalent width, 


W Ic = J (Iq - I(v)) dv 
Line 


(E6) 


Dividing first by and then by Equation E6 becomes 

W = (1“'!>(n))/4i du = 

J 

Line Line 




D(v) dv 


(E7) 


t. 


Frcra Equation E5» 
W s (it>' 

or 

W = (V 

Equations E5 and 
adjusted, once a 
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this is equivalent to 
/ ^ 

0 0 

Line 




D'(v) dv 


/^J) W 
0 0 


(E3) 


£8 allow the measured line intensities to be 

suitable value for umbral intensity ^ is assumed, 
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